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Abstract 


An  Analytical  study  was  conducted  to  determine  the  natural  frequencies  and 
mode  shapes  for  laminated  anisotropic  plates,  including  the  effects  of  shear 
deformation  and  rotatory  inertia,  by  using  the  Galerkin  Technique.  Three  different 
boundary  conditions*  simply-supported,  clamped,  and  two  opposite  sides  clamped, 
two  opposite  sides  simply-supported,  were  considered.  Two  different  graphite- 
epoxy  symmetric  plates  were  used  in  the  analysis.  Convergence  characteristics  and 
the  effects  of  length  to  thickness  ratios  were  investigated.  Comparison  to 
classical  results  and  contour  plots  for  several  mode  shapes  are  provided. 

It  was  found  that  as  the  length  to  thickness  ratios  were  reduced,  shear 
deformation  effects  significantly  lowered  the  natural  frequencies.  Analysis  also 
showed  that  rotatory  inertia  effects  were  very  small.  Convergence  characteristics 
for  all  three  boundary  conditions  were  very  good  and  excellent  agreement  with 
classical  solutions  was  achieved. 


In  recent  years  interest  in  the  use  and  application  of  composite  materials 
has  greatly  increased.  This  is  due  in  part  to  their  high  strength  to  weight  ratios 
and  the  fact  that  they  can  be  tailor  made  for  specific  applications.  Because  of 
their  unique  properties)  they  have  opened  up  numerous  fields  of  research  and 
analysis. 

Past  developments  have  shown  that  the  dynamic  response  of  composite 
plates  departs  more  from  Classical  Thin  Plate  Theory  than  isotropic  ones  do. 
Classical  Thin  Plate  Theory  is  based  on  the  assumption  that  plane  sections  remain 
plane  after  deformation  occurs.  It  has  been  found  that  vibration  analyses  based  on 
this  theory  yield  frequencies  that  are  too  high.  Therefore)  to  gain  better 
•  agreement  with  reality(  the  theories  used  to  analyze  the  response  of  composite 

plates  need  to  include  the  effects  of  shear  deformation  and  rotatory  inertia. 

A  number  of  theories  including  shear  deformation  and  rotatory  inertia  have 
been  proposed  to  date.  Mindlin  Cl  13  introduced  a  two  dimensional  theory  of 
flexural  motion  for  isotropic  elastic  plates.  This  theory  is  based  on  the  fact  that 
the  change  of  displacement  results  from  two  rotations  due  to  bending  and  two 
rotations  due  to  shear  deformation.  This  theory  assumes  that  no  warping  of  the 
plane  section  due  to  shear  occurs>  but  does  include  a  correction  factor  to  account 
for  this  inconsistency. 

Yang>  Norris(  and  StavsKy  C2I3  extended  this  theory  (commonly  referred  to 
as  the  YNS  theory)  to  laminates  consisting  of  an  arbitrary  number  of  bonded 
anisotropic  layers.  They  considered  the  frequency  equations  for  the  propagation 
of  harmonic  waves  in  a  two-  layert  infinite)  isotropic  plate. 
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Whitney  end  Pagano  C 19 3  applied  the  YNS  theory  to  laminated  plates 
consisting  of  an  arbitrary  number  of  bonded  anisotropic  layers)  each  having  one 
plane  of  symmetry  parallel  to  the  central  plane  of  the  plate.  They  employed  this 
theory  to  study  the  cylindrical  bending  of  antisymmetric  cross-ply  and  angle-ply 
plate  strips  under  sinusoidal  loading  and  presented  a  closed  form  solution  for  the 
free  vibrations  of  antisymmetric  angle-ply  plate  strips.  Following  Whitney  and 
Paganoi  Bert  and  Chen  [23  presented  a  closed-form  solution  for  the  free  vibration 
of  simply-supported  rectangular  plates  of  antisymmetric  angle-  ply  laminates. 

Finite  element  analysis  of  laminated  plates  including  transverse  shear 
effects  began  with  Pryor  and  Barker  C 1 3  3.  Their  model  was  based  on  Reissner's 
plate  theory  and  was  applied  to  the  cylindrical  bending  of  a  symmetric  cross-ply 
laminate.  Reddy  C 143  developed  a  simple  finite  element  model  based  on  the  YNS 
theory  and  applied  it  to  the  free  vibration  of  antisymmetrici  angle-ply  plates. 
Reddy  and  others  C 133  also  applied  this  technique  to  orthotropic  laminates  of 
bimodulus  materials. 

Up  to  this  pointi  all  of  the  methods  discussed  could  not  be  applied  to 
laminates  which  possessed  bending/torsional  stiffness  parameters.  This  implies 
coupling  of  the  equations  of  motion.  Recently!  Sathyamoorthy  and  Chia  C 143  used 
the  nonlinear  von  Karman  equations  to  develop  the  theory  to  study  the  large 
amplitude  vibrations  of  anisotropic  skew  plates  for  simply-supported!  damped!  and 
clamped  simply-supported  boundary  conditions.  This  theory  did  include  the 
bending  'torsional  stiffness  parameters.  They  solved  these  equations  using  the 
Galerkin  method  and  the  Runge-Kutta  procedure!  however!  their  applications  were 
limited  to  homogeneous  laminated  plates. 

Thusi  a  need  exists  for  a  method  that  will  take  into  account  shear 
deformation  and  rotatory  inertia  effects  that  includes  the  bending/torsional 
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stiffness  parameters,  yet  it  not  burdened  by  large  computational  requirements. 


Objtctivts 

The  purpose  of  this  thesis  it  threefold.  First,  it  will  prtttnt  a  mtthod  to 
determine  the  natural  frequencies  and  mode  thapet  for  anisotropic  symmetric 
laminated  plates,  including  the  ejects  of  shear  deformation  and  rotatory  inertia, 
using  existing  theories  and  techniques.  This  method  will  include  the  use  of 
bending/torsional  stiffness  parameters.  Second,  the  method  will  be  used  to  analyze 
the  effects  of  shear  deformation  and  rotatory  inertia  on  two  different  symmetric 
laminated  plates.  And  third,  comparison  of  this  study  with  classical  solutions, 
where  available,  will  be  accomplished  to  validate  the  method. 

Approach 

The  approach  used  for  achieving  the  desired  objectives  is  straightforward. 
The  motion  of  an  anisotropic  symmetric  laminated  plate  will  be  modelled  using  the 
YNS  extension  of  the  Hindlin  Plate  Theory  C2#3  and  the  resulting  differential 
equations  of  motion  will  be  solved  using  the  GalerKin  Technique.  To  do  this, 
assumed  functions  will  be  selected  for  each  boundary  condition  and  the  Galerkin 
equations  will  be  established.  From  these  equations  the  eigenvalue  problem  will 
be  formulated.  A  computer  program(s)  will  then  be  written  to  compute  the 
stiffness  and  mass  (or  inertia)  matrix  elements  and  to  solve  the  eigenvalue 
problem.  The  solution  of  this  problem  will  yield  the  desired  natural  frequencies 
and  coefficients  to  determine  the  mode  shapes.  Comparison  of  these  results  with 
classical  solutions,  where  available,  will  be  accomplished.  Also,  to  insure  valid 
results,  a  study  of  the  convergence  characteristics  will  be  performed  by  increasing 
the  number  of  terms  in  each  GalerKin  equation. 


Shear  deformation  and  rotatory  inertia  will  be  studied  by  computing  the 
natural  frequencies  over  a  range  of  length  to  thickness  ratios.  These  frequencies 
will  then  be  compared  with  classical  results.  This  will  first  be  done  without 
rotatory  inertia*  thus  determining  when  shear  deformation  becomes  important. 
Then,  rotatory  inertia  will  be  included  and  the  frequencies  will  be  recalculated 
over  the  same  range  of  length  to  thickness  ratios*  thus  ascertaining  the  effects  of 
rotatory  inertia. 


II. 


Thtory  and  Modelling 

We  will  begin  our  discussion  of  the  theory  used  in  this  thesis  by  starting 
with  some  results  -from  Classical  Laminated  Plate  Theory.  From  there,  we  will 
discuss  the  equations  of  motion  for  a  thick  laminated  plate  using  the  Mindlin  Plate 
Theory,  and  then  discuss  an  approximate  method  of  solution  for  differential 
equations  called  the  Galerkin  Technique.  Galerkin's  technique  will  then  be  applied 
to  our  differential  equations  of  motion  for  three  different  boundary  conditions. 

Anisotropic  Thick  Plate  Theory  with  Rotatory  Inertia 

Classical  Laminated  Plate  Theory  (commonly  referred  to  as  CLPT) 
incorporates  constitutive  relationships  for  an  orthotropic  lamina  through  the  plate 
thickness  resulting  in  expressions  which  approximate  force  resulants  in  terms  of 
displacement  functions.  The  concepts  from  CLPT  are  essential  for  the  later 
development  of  the  equations  of  motion.  We  will  begin  by  describing  the  basic 
constitutive  relationships  for  an  individual  lamina.  The  reader  should  refer  to 
References  C13  and  C83  for  an  in-depth  development  of  these  relations. 

The  basic  constitutive  relationships  for  a  single  orthotropic  layer  in  the 
fiber  oriented  reference  system,  as  shown  in  Figure  2.1,  are: 


Figure  2.1  Definition  of  Coordinate  System 
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where  €  j,  €  2»  ind  f 3  *r*  th*  normsl  strAins,  <5»  *nd  Art  th*  sh*Ar  str*ins, 

o  . ,  O-i  And  O-  Art  th*  normsl  stress**!  And  o  o_ .  And  O.  Art  th*  thttr  *tr*t**». 
12  3  4  5  6 

The  S.j  term*  form  th*  compliAnc*  mttrix  And*  *xpr*t**d  in  term*  of  th* 
engineering  con*t*nt*,  a r*: 
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where  8.  ere  Young's  moduli  in  the  ith  direction!  V..  i*  Poisson's  rstio  for 


tr*nsv*rse  strain  in  th*  jth  direction  whan  stressed  in  the  ith  direction  And  G.^  it 
the  shear  modulus  in  th*  i-j  pl*n*. 


If  w*  invert  Bq.  (1)  to  g*t  stresses  in  terms  of  strains  w*  h*ve 


<0>  ■  CQ'3  (€> 


where  CQ'3  is  th*  reduced  stiffness  matrix  end  hts  the  form 
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If  we  wish  to  rotate  our  principal  lamina  axis  an  amount  d  (in  the  1-2  plans) 
with  rsspsct  to  soms  geometrical  axis  (refer  to  Figure  2.1)  then  our  stiffness 
matrix  must  under  go  a  tensor  transformation.  The  transformation  matrix  would  be 
defined  as 
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where  m»cos(6)  and  n«sin($).  The  transformed  stiffness  matrix  would  then  be 


CQ'J  «  CT3  CQ'3  mT 
and 


(7) 


(O)  ■  CQ'3  «> 


(8) 


At  this  point  we  shall  assume  for  a  thin  laminae  a  state  of  plane  stress. 
Therefore*  O  ■  4,  From  Bq.  (8)»  we  have 


<z  •  <Q'i3/QV<x  4  (Q'23/Q,33Ky  *  «W«y 


Substituting  Bq.  (11)  back  into  Eq.  (8)  the  resulting  CQ'3  matrix  becomes  a  5x5 
matrix  referred  to  as  CQ3. 
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We  can,  therefor*!  writ*  Bq.  <8)  ti 
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We  now  with  to  build  a  laminate  from  N  perfectly  bonded  lamina  and  exprett 


the  forces  and  moments  acting  on  this  laminate  in  terms  of  displacement  functions. 
Before  we  do  this»  we  need  to  make  an  assumption  concerning  for  the  laminate 
and  we  need  to  discuss  our  sign  convention. 


We  *h»ll  assume  that  This  implies  that  a  lint  perpendicular  to  the 

midplane  will  not  stretch  under  deformation  and  is  an  accepted  inconsistency  in 
plate  theory.  In  reality*  €  is  not  zero*  but  is  small  compared  to  the  other  strains. 
For  our  laminate*  it  means  there  will  be  discontinuities  in  at  the  lamina 
boundaries  (they  too*  will  be  small).  With  this  assumption  we  can  assume  a 
displacement  field  of  the  form 

u  *  u°(x,y,t)  +  z»x(x*y*t) 

v  «  v°(x,y,t)  +  z?y(x,y,t)  (17) 

w  ■  w(x*y,t) 


where  u,v*  and  w  are  the  x,y,  and  z  coordinate  displacements  respectively*  u°  and 

v°  are  the  displacements  of  the  midplane  of  the  laminate,  and  *  and  ▼  are 

x  y 

rotations  of  a  line  perpendicular  to  the  midplane  due  to  bending.  It  is  important  to 

note  here  that  in  this  thesis  f  and  f  will  be  thought  of  as  rotations.  Thus,  for 

x  y 

an  axis  system  as  defined  in  Figure  2.2,  a  right  hand  rotation  about  the  positive  y 

axis  will  give  us  a  positive  *  rotation.  On  the  other  hand,  we  will  think  of  w 

x  *x 

and  w  (the  commas  denote  differentiation)  as  rates  of  change  of  displacement 
with  respect  to  the  appropriate  coordinate.  Some  authors  look  at  the  rotations  v 

x 

and  ?  as  rates  of  change  of  displacements  also  (see  Ref.  C133).  Depending  upon 
the  point  of  view  one  takes*  the  sign  on  these  rotations  is  different.  However,  if 
one  is  consistent  throughout  his  development  with  his  notation  then  either 
approach  will  yield  the  same  result  related  to  generalized  functions  such  as 
eigenfunctions  and  eigenvalues. 


At  this  point  we  have  not  modelled  sheer  deformation  and  one  should  take  note 
that  the  above  curvatures  do  not  include  their  effects.  Now  Bq.  (15)  can  be  written 
(for  the  kth  layer  of  a  laminate)  as 


(23) 


The  resulant  forces  and  moments  acting  on  a  laminate  are  obtained  by 
integrating  the  stresses  in  each  layer  through  the  laminate  thickness.  Thus,  we 
have  (see  Figure  2.3  for  the  geometry  of  an  N  layered  laminate) 


Figure  2.3  Geometry  of  an  N  Layered  Laminate 
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Tht  integrations  in  Sqs.  (24)  And  (25)  can  be  rearranged  to  take  AdvAntAg*  of 

the  fact  thAt  the  stiffness  (nAtrix  tor  a  lamina  is  constant  within  the  lamina.  We 

CAn  Also  recAll  thAt  €°  ,  <°  ,  Y°  ,  K  t  K  »  And  ic  Are  not  functions  of  2  but  are 

x  y  xy  x  y  xy 

middle  surface  velues  of  the  lsminAte  And  can  be  removed  from  under  the 
summation  signs.  Therefore,  Bqs.  (24)  and  (25)  can  be  written  as 


(26: 


<  28  J 


Aij  “  <Zk"V,) 


siJ  *'1j1<5iA  <2«2-iic-iJ> 


(2?) 


°ij  *  1/3£<5U>i>  <2k3-!k-l3> 


(30) 


In  Bqs.  (28),  (29),  and  (30)  the  A. ere  called  the  extensional  stiffnesses,  thf 

B.  's  art  called  the  coupling  stiffnesses,  and  the  D.'s  are  called  the  bending 
ij  l  j 

stiffnesses. 


We  now  need  to  determine  an  expression  for  the  transverse  shear  forces  Q 
and  Qy  in  terms  of  displacement  functions.  To  do  this  we  need  to  discuss  Hindlin 
Plate  Theory. 

In  classical  plate  theory  the  Kirchhoff  hypothesis  states  that  straight  lines 

perpendicular  to  the  neutral  axis  in  the  undeformed  state  remain  straight  and 

perpendicular  to  the  neutral  axis  after  deformation.  Hindlin  discarded  this 

assumption  and  assumed  that  these  lines  would  remain  straight  (no  warping)  but 

would  not  remain  perpendicular  to  the  midplane  after  deformation.  The  assumption 

of  no  warping  is  not  correct  and  he  did  introduce  a  factor,  k,  as  a  means  of 

compensating  for  it.  Therefore,  the  slope  of  the  midplane  (w  and  w  )  now 

»x  »y 

consists  of  a  rotation  due  to  bending  (<r  ,  e  )  and  a  rotation  due  to  shear 

x  y 

deformation  0^,  V  ).  See  Figure  2.4. 
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From  the  definition  of  engineering  strein  we  hive 


Y  •  u  ♦  w 

XZ  »Z  ,x 


Substituting  the  first  expression  in  Bq.  (1?)  into  Eq.  (31),  we  hive 


Therefore 


U#i  "  Tx 


Y  •  w  ♦  f 

XZ  ,X  X 


ind  from  similier  reisoning 


Y  ■  w  +  »  (34) 

y*  »y  y 

Here  we  cm  see  thit  the  curviture  is  i  function  of  both  bending  ind  sheer 
deformetion.  Substituting  Eqs.  (33)  end  (34)  into  Bq.  (16)  end  introducing  k,  the 
correction  fictor  for  issuming  the  sheer  werping  to  be  e  streight  line*  we  hive 


°44  °45  {V'y 


*45  °55  )W,x+,x 


The  sheer  forcesf  Q  end  Q  mey  now  be  determined  by  integnting  the  sheer 
stresses  end  y  through  the  limine  thickness  end  summing  over  the  leminete. 

We  hive 


Q  »  T 

*  •'-h/2  xz 
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As  before  the  integration  can  be  rearranged  to  take  advantage  of  the  fact 
that  the  stiffness  matrix  for  a  laminate  is  constant  within  the  lamina.  We  can  also 


note  that  w  ,  w  ,  *  ,  and  f  are  not  functions  of  2  but  are  midplane  values  of  the 
,x  ,y  x  y 

laminate.  Therefore*  we  have 


Com 

y 
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A44  A45 

h 

A45  A55_ 

fw  +T  *' 

,y  y 

]  w  +<r 

L  »x  x) 


(38) 


where  the  A.,  terms  are  defined  in  Eq.  (28).  We  now  have  the  required  forces  and 
U 

moments  in  terms  of  displacement  functions  as  we  desired  and  may  proceed  on  to 
the  governing  equations  of  motion  for  a  plate. 


The  governing  equations  of  motion  for  a  plate  may  be  derived  by  formulating 
the  Lagrangian  function  (in  terms  of  plate  variables)  for  the  plate  and  applying 
Hamilton's  principle  to  that  Lagrangian.  This  approach  will  also  yield  the  required 
boundary  conditions.  Reference  T63  formulates  the  Lagrangian  for  a  plate 
including  the  effects  of  transverse  shear  and  rotatory  inertia  (modelled  using  the 
Mindlin  Plate  Theory).  Following  the  derivation  from  that  text,  we  have,  after 
applying  Hamilton's  principle  and  Green's  theorem,  the  following: 
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(39) 


where  the  double  integral  over  the  domain  represents  the  equations  of  motion  for 
the  plate  and  the  line  integral  around  the  contour  represents  the  required 
boundary  conditions.  Also 


To  determine  the  equations  of  motion  tor  the  plate  at  any  tim*  tf  wt  take 
the  double  integral  expression  over  the  domain  and  note  that  the  variation  cannot 
be  zero  over  this  region.  Therefore*  the  coefficients  for  each  variation  coordinate 
must  be  equal  to  zero.  Thus  we  have 


M 

♦  M  -  Q 

It 

x  ,x 

xy,y  x 

X 

M 

♦  M  -  Q 

I? 

xy,x 

y,y  y 

y 

Q 

♦  Q  „  +  q  ■ 

Pw 

x,x 

yiy  ^ 

If  we  assume  the  time  dependence  to  be  harmonic)  then  we  can  separate  out 
the  time  variable  in  Bq.  (41)  and  we  have 


M  ♦  M  -  Q  ♦  -  e 

x,x  xy,y  x  x 


<42a> 


M  ♦  M  -  Q  +  u  It  *  0 

x/,x  y,y  y  y 


<42b> 


Q  *  Q  ♦  N  w  +2N  w  ♦  N  w  ♦  pu  w  *  0  (42c) 

x,x  y,y  x  ,xx  xy  ,xy  y  ,yy 

where  0)  is  the  frequency  of  vibration  and  we  have  let  q  •  N  w  ♦  2N  w  * 

^  x  »xx  xy  fxy 

Nyw  .  This  will  allow  the  use  of  our  equations  to  solve  the  linear  bifurcation 
problem.  We  will  retain  this  expression  throughout  our  equation  development  but 
will  not  consider  it  when  solving  the  vibration  problem.  Thus#  Bqs.  <42a)  thru  (42c) 


are  our  equations  of  motion. 


Before  we  proceed  we  must  determine  the  appropriate  boundary  conditions 
for  the  system  of  equations  which  are  sufficient  to  assure  a  unique  solution. 


These  can  be  determined  from  the  line  integral  expression  in  Eq.  (39).  The 
boundary  conditions  associated  with  Eq.  (42a)  are  represented  by  and  its 
coefficient.  Thus,  we  have 


dy  +  M  dx)  &¥  *  8 

x  '  xy  x 


(43a) 


Similiarly,  the  associated  boundary  conditions  for  Eqs.  (42b)  and  (42c)  are, 
respectively 


fr 

Sr" 


<M  dx  -  M  dy>  • 

y  xy  y 


Q.  dy  +  Q  dx)  =  • 
x  y 


(43b) 


(43c) 


We  now  have  the  equations  of  motion  and  the  required  boundary  conditions 
and  can  proceed  with  finding  a  solution  for  this  system.  Before  we  do  this, 
however,  we  need  to  apply  a  restriction.  This  thesis  will  only  be  concerned  with 
symmetric  laminates.  Therefore,  all  of  the  coupling  stiffnesses  (Bj/s)  are  equal 
to  zero  and  the  extensional  stiffness  term,  A^g,  is  also  equal  to  zero.  With  this 
restriction  Bqs.  (2?)  and  (38)  become 
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Sine*  we  have  Arrived  at  a  set  of  partial  differential  equations  that 
describe  the  motion  of  our  plate*  we  now  need  to  solve  them.  This  thesis  will  use 
An  approximate  method  celled  the  Galerkin  technique.  As  opposed  to  a  technique 
like  the  Ritz  method*  this  technique  involves  the  direct  use  of  the  differential 
equetion  end  does  not  require  the  existence  of  a  functioned  Thus*  it  hes  a  broader 
range  of  epplicetion  then  the  Ritz  method.  However*  in  the  eree  of  solid 
mechenics*  the  two  ere  closely  releted. 

The  bssic  idee  behind  the  Galerkin  technique  cen  be  discussed  quite  briefly. 
Suppose  we  went  to  solve  the  equetion 

L(u)  =  0  <4 A) 

where  L  is  some  differentiel  operetor  in  two  veriebles  whose  solution  setisfies 
homogeneous  boundery  conditions.  We  will  look  for  an  Approximate  solution  in  the 
form 

N 

u(x ,y)  -  I  c .  <x ,y)  <4?) 

i-1  1  ' 

where  f^ix.y)  is  a  system  of  functions  which  satisfy  the  boundary  conditions  and  c. 
are  undetermined  coefficients.  We  will  assume  the  functions  4^(x*y)  to  be  linearly 
independent.  It  is  important  to  note  here  that  if  we  are  to  get  an  "accurate” 
answer*  the  system  of  functions  ♦  .<x»y>,  must  be  complete  in  the  given  region.  If 
they  are  not*  then  we  are  excluding  part  of  the  solution  and  our  results  will  be 
very  misleading.  If  u(x*y)  is  to  be  an  exact  solution*  then  L(u)  will  be  identically 
equal  to  zero.  If  LCu)  is  continuous  then  we  are  saying  that  L(u)  is  orthogonal  to 
all  of  the  functions  of  the  system  ♦  .(x*y).  Mathematically  speaking  we  have 


f  f  N 

JD  L<u(x,y>)#.<x,y>  dxdy  -  JD  L<  2  c  .  t  .  <x  ,y))4  (x  ,y)  dxdy  *  0  <48> 

1  j*l  J  J 

where  i  •  l„.„N  individual  equations.  From  •this  system  of  equations  we  can  solve 
for  the  undetermined  coefficients  and  therefore  arrive  at  our  solution. 

The  question  now  arises  what  happens  if  our  assumed  function  does  not 
satisfy  the  boundary  conditions.  If  it  doesn't)  we  want  to  force  our  approximate 
functions  ♦  ,(xly)  to  satisfy  them  also.  Thus«  for  the  case  of  two  variables  we 
would  set  up  a  GalerKin  equation  as  a  line  integral  around  the  boundary  that  would 
force  the  function  to  satisfy  the  boundary.  For  homogeneous  boundary  conditions 
this  line  integral  in  and  of  itself  would  also  equal  zero. 


23 


At  this  point  we  art  ready  to  formulate  the  Galerkin  aquations  tor  the  plate 
equations  ot  motion.  Our  approach  will  be  to  formulate  the  Galerkin  equations! 


make  appropriate  substitutions!  and  then  normalize  the  equations. 

In  formulating  the  Galerkin  equations  we  must  be  careful.  If  we  could  find 
approximate  functions  which  would  satisfy  both  the  geometric  and  natural  boundary 
conditions!  then  we  could  apply  Sq.  (48)  to  each  of  our  three  equations  of  motion 
and  our  Galerkin  equations  would  be  formulated.  Howevert  finding  such  functions 
is  extremely  difficult!  if  not  impossible.  We  canf  however!  find  functions  which 
satisfy  the  geometric  boundary  conditions.  This  does  impose  the  added  condition 
that  we  now  must  formulate  a  Galerkin  equation  for  the  boundaries  as  well  as  the 
equations  of  motion.  We'll  start  by  formulating  the  Galerkin  equations  for  Sqs. 
(42a)  and  (43a).  If  we  assume  that  is  the  approximate  rotation  function  about 
the  y  axis  and  use  E!q.  (48) .  we  have  for  Eq.  (42a) 


a 


(M  ♦  M  -  Q  +  «2I?  )  9'  dxdy 
x,x  xy,y  x  x  x 


(49) 


where  is  made  up  of  a  series  of  functions  (each  which  satisfy  the  geometric 

boundary  conditions)  where  each  term  in  the  series  is  multiplied  by  an 

undetermined  coefficient!  A  .  T'  on  the  other  hand!  it  a  single  term  of  the  T 

mn  x  x 

approximate  series  and  has  no  undetermined  coefficient  associated  with  it.  This  is 
a  result  of  the  Galerkin  derivation  and  several  examples  of  this  may  be  found  in 
References  CS3  and  C173.  For  each  term  in  the  approximate  ?  series!  a  new 
equation  will  be  generated.  We  should  note  here  that  one  of  the  conditions  for 
solving  for  the  undetermined  coefficients  is  that  there  are  as  many  equations  as 
there  are  undetermined  coefficients.  Thus!  if  we  have  men  undetermined 
coefficients  in  our  approximate  seriesi  T  .  Bq.  (49)  will  generate  men  equations. 


We  will  now  turn  our  attention  to  the  formulation  of  the  Galerkin  equation 


for  the  boundary,  Eq.  (43a).  Eq.  (43a)  if  different  from  Eq.  (42a)  in  that  it  still 
contains  the  variation  of  ?x*  This  poses  no  real  problem  in  that  the  Galerkin 
method  can  be  related  to  the  functional  expression  in  Eq.  (39).  Using  Eq.  (43a)  and 
following  the  procedure  from  Reference  [5  3  we  have 

J(-M  dy  ♦  M  dx>  ?'  *  0  (50a) 

r  x  «y  x 

where  ?'  is  again  a  single  term  of  the  ?x  approximate  series  and  does  not  contain 
an  undetermined  coefficient.  The  undetermined  coefficients  are  hidden  in  the 
and  terms  (see  Eq.  (44)  where  is  replaced  with  V^).  If  we  use  the  plate 
coordinates  (see  Figure  2.2)  to  define  the  contour  for  the  integral  in  Eq.  (50a),  it 
becomes  (see  Ref.  C183) 


<M  <0 ,y)t 
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x<*,X> 


-  Mx<a,y)f'x<a,y>)  dy 


♦  f  <M  <x,0)9'  <x,0) 

J0  xy 


-  M  <x ,b)f '  (x ,b) >  dx 
xy  ’  x 


0 


(50b) 


Bqs.  (49)  and  (50b)  form  the  Galerkin  equations  for  the  first  coordinate  <r  .  In 
order  to  reduce  the  size  of  the  eigenvalue  problem  we  will  be  solving  later,  we  can 
approximate  the  two  Eqs.  (49)  and  (50b)  by  combining  them  into  one  equation.  We 
have  therefore 


<h  ♦  M  -  Q  ♦  u2If  }  V  dxdy 
x ,x  xy ,y  x  x  x 

(Mx<0,y)  ?'(0,y)  -  Mx<a,y)  ?'<a ,y)3  dy 
<M  <x,0)  f'(x,0)  -  M  <x,b)  ?'  <x  ,b) )  dx 


0 


(51) 


Now  substituting  Bq.<44)  *nd  (45)  into  (5i)  end  combining  terms  we  heve 
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We  choose  to  normelize  Sq.  <52>  using  the  normelizetion  scheme  simiiier  to 
the  one  used  in  Reference  C183. 
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(53) 


where  d^  end  ere  the  normelized  bending  end  extensionel  stiffnesses 


respectively*  is  Young's  Modulus  in  the  direction  normel  to  the  lemine  fibers*  R 
is  the  espect  retio*  s  is  the  length  to  thickness  retio*  to  is  the  normelized 


frequency,  p  it  the  product  of  lamina  density  and  lamina  thickness  tunned  over  th* 
laminate*  %  it  the  normalized  x  coordinate*  "’l  it  the  normalized  y  coordinate*  and  h 
it  the  plate  thicknttt.  We  will  alto  make  the  aubatitution 


-h/2 

[  n  dz 

J  -h/2 


where 


Ph3/12  »  ph2/12 


N 


(54) 


p  ■  2  P,etk 

k-1  *  K 


(55) 


The  reaton  for  this  aubstitution  is  that  it  will  allow  ut  to  build  up  a  laminate  o-f 
different  materials.  Thus*  we  can  use  our  homogeneous  theory  to  handle  hybrid 
composite  plates.  After  these  substitutions  and  subsequent  simplifications*  Bq. 
(52)  becomes 
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Wt  are  now  ready  to  formulate  Galtrkin's  equation  for  Eq».  (42b)  and  (43b). 
Using  tht  same  procedure  as  before#  we  arrive  at 


<M  ♦  M  -  Q  +  W2If  }  ?'  dxdy 
xy,x  y,y  y  y  y 


I*  D 

♦  J  <Mxy(0,y)  iy<0,y)  -  Mxy<*,y>  *'<»./>>  dy 


J(M  <x,0)  a'fx.S)  -  H  <x,b)  ^'(x.b)} 
e  *  y  y  y 


dx  ■  0 


where  <r  it  tha  assumed  rotation  function  about  tha  x  axis.  At  before  tha  doubla 

y 

intagral  will  force  our  approximata  rotation  function  to  satisfy  tha  plata  domain 
whila  tha  two  lina  integrals  will  force  the  function  to  satisfy  the  natural  boundary 
conditions.  Now  substituting  Bqs.  (44)  and  (45)  into  (57)  and  combining  terms  wa 


f‘fb  <DU?,,x,  *  <°124'V\,XX  *  °2i\,ry  4  °M?y,x,  4  2t,2i"y,yy 
4  °22’y,yy  ‘  “mV  ‘  Kfl«5,y  *  V  K  **  dy 

*X  <0l4?x, ,<*•>'>  *  ■>24’ytx(''y>  4  <\.y"'y,4,y, ,">•>'»>  *;«.>-» 


-  10ll’x, ,U'y>  *  D24V, y<*’y>  4  E>.4<%,y<4,y>44y,x<4,y)' 5  V*'y>  dy 


*/*  ",l2’x,x<x',)  4  0 


<x,0)  ♦  D,,<t  <x,0)+a  (x,0))}  f  <x,0> 

22  y,y  ’  24  *,y  '  y,x  1  y  ' 


"<Dl2fx,x<x,b>*D22Tyly<x,b>4024<,x  y<x»bH%  x<x,b))}ay(x  ,b)dx  -  0  (58) 


If  we  normalize  Sq.(58)  using  the  same  normalization  scheme  as  before*  we 


fX  “ 


16%, «  +  <d12 +  d66)%*44  *  R2d2 6%,W  +  d66%,4!;  +  2Rd26%(*i) 


♦  ^22%, m  "  ks2a44%  '  ksa44W,1  +  <“2/12)V  d%  dD 

"J1  ldU,%.%<l»,>  +  Rd26%,V0>1))  +  %<0’1” 

^  0 


<duV%<l,,)  +  Rd24?1,V1»,>  +  tWR%,V1'1|)*?M<l,1'>>5  d1* 


♦  Rd225iJ ,n(%,0)  +  d26<R%,’)<*,0>+%I4<*,0))>  V*’0> 


"{d12’^l^<%,1)4Rd22’tlV%'1) 


+d26<R%,,»<',1)+%,4<',1>)>V*,1)  d*  "  0 


We  now  need  to  formulate  Galerkin's  equation  for  the  last  equations  of  our 
system*  Bqs.  (42c)  and  (43c).  We  have 


♦  Q  ♦  N  w  +2Nw  +Nw  +  Pt#2w)  Ul'dxdy 


x,x  y,y  x  ,xx  xy  ,xy  y  ,yy 


/• b  _/ 

<Q  <*,y)  U<i ,y)  -  Q  (a,y)  U(a,y)>  dy 

J  0  * 

/•*  _/ 

♦  (Q  <x  ,8)  U(x,0)  -  Q  <x,b)  W<x,b)>  dx  «  0 

J  8  y  Y 


where  W  is  the  assumed  displacement  function.  Again*  the  double  integral  is  to 
satisfy  the  domain  and  the  two  line  integrals  are  to  satisfy  the  boundary. 
Substituting  Bqs.  (44)  and  (45)  into  (64)  and  combining  terms  we  have 


2? 


f‘f*  <kA55’x,x  *  kfl55“,xx  *  M^y.v  *  kA«“,xx  '  k2No“,xx 
*  2k/*o“,xx  ‘  klNo“,yy  *  w2u)  dxd>' 


fb  - 

<kA55?x<0,y>  ♦  kA5^w>x(8,y)>  W<0,y> 

w  0 


-  <kA55fx<a,y)  ♦  kA^w  x<a,y>>  W<a,y>  dy 

/•*  _/ 

+1  <kA..T  (x,0)  +  kA..w  <x,0»  W<x,0> 

Je  r  **  >y 


-{kA-.t  <x ,b)+KA- -w  <x,b)>  W<x,b)  dx  -  0  <61 ) 

y  i  y 

whir#  Ny«-k1No,  Nx— k2NQ,  and  N^-kgN^ 

W*  will  now  normalize  Eq.(61)»  however,  we  need  to  define  two  more 
normalization  factor*.  Thay  ara 


x.  -  N  b2/E_h3 

10  4 


x2  ■  pb4u2/82h3 


Therefore,  Eq.(6i)  will  bacoma 


U4 


„  *  ka„w  ♦  kRsa..?-  _  +  kR  a..w  k  -  k„x.<R  /*  )w 

0  J0  55  55  44  ’l, fl  44  ,%%  2  1  , 

+  2k3<R3/s2)XjW  ^  -  k1(R4/s2)X1w  ^  +  (R4/42)X2w)  w'd^dH 


f 1  _  -/ 

♦  {ksa55f%(0,H>  +  ka55w  4<0fi»)  W<8 , 1 

W  0  * 


-  {ksa55f^(lf^)  ♦  ka^w  W<1  ,1) 


♦J  (ksa44t ,,(%,*)  ♦  Rka44w|tl<%,0))  W<*,0> 

w  0 


-  <K**44»1|<%|1)  ♦  Rka44w  T1<*» 1  >>  uU.l)  di  -  0  <63> 

At  this  point  wi  art  ready  to  pick  a  boundary  condition!  find  approximate 

functions  for  T  ,  »  ,  and  w  for  that  condition!  and  substitute  into  Bqs.  (56),  (59), 
x  y 


and  (63). 


Since  the  Galerkin  equations  have  been  formulated  it  is  possible  to  choose 


the  displacement  and  rotation  functions  (w*  f  *  e  )  which  satisify  the  three 

x  y 

equations  for  each  boundary  condition  being  considered.  These  differential 
equations  can  be  integrated  resulting  in  three  algebraic  equations  (for  each 
boundary  condition).  The  eigenvalue  problem  can  be  formulated  from  which  the 
natural  frequencies  and  mode  shapes  can  be  obtained.  Let  us  begin  with  the 
simply-supported  boundary  condition. 

For  a  plate  simply  supported  on  ill  sides,  we  know  from  plate  theory  that: 


at  x  *  t,i 


w  *  f  * 

y 


(64) 


M  «D..*  ♦D.-f  ♦  Di.(*  +¥  )*• 

x  11  x,x  12  y*y  16  yfx  x,y 


and 


at  y  ■  6  ib 


w  ■  T  ■  • 
x 


(65) 


My  '  D12  *  D22  *  D2i  ‘V  *  ’»„>  *  ‘ 


Therefore!  we  shall  choose  for  our  admissible  functions  (functions  which  satisfy 
the  geometric  boundary  conditions  and  are  continuously  differentiable  one  time) 


e 

x 


U  - 


I  X  A_ 


ra*l  n*l 


ran 


X  X 


m=l  n=*l 


ran 


1  1  Cmn 

s-1  n-1  1,0 


cos(mfx/a)  sin<n<y/b) 


sin(ratx/a)  cos<nty/b> 


sin(ra<x/a)  sin<n<y/b> 


(66) 


or  in  terms  of  normalized  coordinates 


».  *  Z  Z  A  cos<mt^)  sin<nt‘!l) 
a  .  ,  mn 

m=>l  n»l 


?_  »  Z  ZB  sin(mt%)  cos(nt1l) 
’l  .  ,  mn 

ro=l  n«l 


<6? ) 


W  »  Z  Z  C  sin(mft)  sin(ntt) 

.  _ .  mn 

m*l  n"l 

In  order  to  substitute  into  Eqs.  (56)*  (59)»  and  (63)  the  following  derivatives 
are  needed  (note:  all  ZX  are  over  m»l  to  *»  and  n»i  to  «►). 


= 

z 

z 

-A 

mn 

at 

z 

z 

-A 

mn 

a 

z 

z 

-A 

mn 

Vi 

a 

z 

z 

Amn 

Vll 

* 

z 

z 

An 

v% 

3 

z 

z 

B 

mn 

V** 

z 

z 

mn 

V<1 

a 

z 

z 

-B 

mn 

Vi 

3 

z 

z 

-B_ 

on 

Vll 

“ 

z 

z 

®mn 

w,% 

■ 

z 

z 

Cren 

w,« 

a 

z 

z 

-C 

mn 

nr 

m 

13 

- 

z 

z 

C 

mn 

".1 

a 

z 

z 

C 

mn 

W,H 

= 

z 

z 

-C 

mn 

ml  sin(nrH)  sin(nt1)) 

m2t2  cos(mtl)  sintnt1)) 
2 

mnt  sin(mt%>  cos(nfn) 

nt  cos(mtt)  cos(nfn) 

n2t2  cos(i»«)  sin(nlH) 

ret  cos(mtt)  costnfl) 

ra2t2  sin (ml|)  cosCnt1)) 
2 

mnt  cos(mtt)  sin(nH) 

nt  sin(mtt)  sintnt1)) 

n2t2  sin(mt|>  cos(ntt) 

mt  cos<ratt)  sin(nti)) 

m2t2  sin(retfc)  sintnt^) 
2 

rent  cos(mt£)  cos(nttl) 
nt  sin(mt%)  cosCnt11)) 
n2t2  sin(mt$>  5in<ntH) 


<68) 


Thus  from  Bq.  (36) 


•  -  *1-1  2  2 

I  Z  {-d.,A  ra  t  cos<mt§)sin<nt*) 

*1  n-i^eJs  11  mn 


2Rd.  ,A  rant  sin<mt%)co&(nt^> 
1 6  ran 


-  R2d,.A  n2t2cos<mt%)sin(nt’l)  -  d.  .B  ra2t2*i n<nt*)cos<nt*l) 

66  mn  16  mn 

-  R(d, _+d, ,)B  mnt2cos<mt*)sin<nt*>  -  R2d0.B  n2t2sin(rat*>cos<nti» 

12  66  ran  26  mn 

2 

-  ks  a5gAmncos<mtt)*in(nt'Tl)  -  ksiggC^mtcc^mttlsi n<nt'tl) 

+  <52/12)A  cos(rat$)si  nCnt1)) )  (cos<pt%)sin<qt11)  }  d$  d*l 


llAran 


ratsin<0)sin<nftf)  -  Rd 


.  _B  ntsin<B)*in(nt‘,l) 
1 2  ran 


♦  d.  .[RA  ntcos<8)co3<ntTl)  ♦  B  mtco*<0>co»<nt‘n) 3>  {cos(0>»in(qt1))> 

1 6  mn  ran 

-  <-d. .A  ratsin(mt)sin<nt^>  -  Rd. 0B  ntsinCmOsintnt11)) 

1 1  ran  1 2  ran 

♦  d.  .[RA  ntcos<ret)cos<nt1l)  +  B  ratco»<mt)co»<nt1l) ) } 

1 6  mn  ran 

♦  {cos(pt)sin(qt‘,l)>  d1*! 


<-d, .A  mtsi n(mt$) si n<0) 
1 6  ran 


-  Rd^.B..  ntsin<rat%)sin<0) 
26  ran 


♦  d,.[RA  ntcos(mt$)co*<0)  +  B  ratcos(ratl)cos<0) ] }  (cos(ptt)si n(0) ) 

66  mn  mn 

-  (-d. .A  ratsin<ret%)sin<nt)  -  Rd0,B  ntsin<mtl)*in(nt) 

16  mn  26  ran 

♦  d , . [ RA  nlcos(nrft)cos<nt)  ♦  B_  ratcos<rott)cos<nt) ] ) 

66  mn  mn 


•  <co*<pt%)sin(qt>)  d%  ■  0 


<69) 


From  Sq.  (59) 


Z  Z  r»2l2cos(m«%)*in<ntn)-R<d12+d^)Aronmn<2*in<ml^)co*(ntl) 

w*ln»W0w0  mn 

-  R2d26Amnn2t2co*<nrt%)*i"<n«)  -  d66Bmnn>2t2*in(r«t%)co5(n<^) 

-  2Rd~ .B  r»n<2co*<nrU)*in<n<’))  -  R2d_B  n2t2sin(m<%>cos<nt^> 

26  mn  ran 

-  ks2a. .B  si  n<mt%)cos<nt1l)  -  kRsa  C  ntsin(mt%)cos<nt^) 

44  mn  *»**  ran 

+  <u2/12)B  sin(ml%)co*<nt*)>  <*in(p<Ocos<qt1!»}  d%  d1) 
mn 


Jm  i 

{-d1<5Amnmtsin<0>sin<n<1>  -  Rd26B[innlsin(0)sin<nlt^) 

0 


♦  d  [RA  n<cos(0)cos(nt^l)  *  B  nrtco*<0>cos(n<il) ] )  (sin(0)cos(qt1|))> 

66  mn  mn 

-  <-d,  ,A  nrtCsin<nrtOsin(nl^l)  -  Rd9XBmnnlsi n<mt)si n(n<i)) 

I  o  mn 

♦  d,,[ftt  nfcos(raf )cos<nt*l)  ♦  B  mtco*<mt)cos<n«>]} 

66  mn  mn 


*  {sin(p<)cos<q<1H))  d1) 


;:<-v 


A  mt*in<mH)sin<0)  -  Rd,_Bitvnn<sin(m<%)sin<0) 
!  mn  22  mn 


♦  d„.  tfto  n*cos<rel<>cos<8)  +  B„  nri(co*<nrt(%>co»<0>3>  {sin<p<|)cos(0>> 

26  mn  mn 

-  <-d12AmnnriCsin<mtt)sin<nl)  -  Rd22B|nnnt*in<mH)sin<n1() 

♦  d„,IIW  n<cos(m**>cos<nf>  +  m<co*<m*0>co*<n*>  J) 

26  mn  mn 


«  {»in<p<%)cos<ql)}  <5%  ■  0 
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I 


And  fr om  Sq.  (63) 


Z  iff  (-k*a__A  „nrt(sin(mt%)sin<nO)  -  ka„C  m2<2sin(m<*)sin<nt4) 
r.1  n-1 JbJb  55^ran  55  mn 


-  kRsa-.B  ntsin<»t^)5in<nt1))  -  KR^i„*C  n^t^*in(nrt^)sin(n<1)) 
44  nn  44  mn 

♦  k_x. <R2/s2)C  i»2t2sin<mt*)sin<n<n) 

£  i  mn 

♦  2k0(R3/*2)X, C  mnt2co*<nt%)cos(nH> 

o  l  mn 

+  k  <R4/s2)X.C  n2l2sin(ml$)sin<nti)) 
i  l  mn 


*  <R  /»‘c)x2Cmn*in(m<(l)ain(nt1))}  <sin<pl%)sin<qt^)>  d%  d1) 


/.  ‘"”5^ 


icos<0)si  n(nt1))  ♦  kaggCranmtcos(0)sin(n<1)))  {*.in(0)»in(q<11))  ) 


-  tksa,-„A  cos(mt)*i  nCnfl)  ♦  ka__C  m<cos<mt)sin(nt1))> 
55  nm  55  mn 


*  (*in<pt)*in<qt1l)}  d1) 


+  I  t*<sa..B  »in<mlli)co»(0)  •»  Rka.-C  nKsi  n(ra<  |)cos<0) )  (si  n(pt%)si  n<0)  } 

J g  44  mn  44  mn 


-  <ksa . .B^*in(mHl)co*(n<)  ♦  Rka„„C  nisi n<mt$)cos(n<) ) 
44  mn  44  mn 


*  <sin<pl^)sin(qt)}  di  ■  0 


The  indicated  multiplication*  are  performed  after  noting  that  sin(0)  * 
*in(p<)  «  sin(qK)  *  •  (p  and  q  are  integer*)  and  co»(0)  »  i. 


>VV- 

s 

V.S"  N 

^v: 
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From  Eq.  (69) 


"  *  /**(**  2  2 

I  Z  I  I  -d.^A  m  t  costmtq>cos(p<$)sintnt1l)sin(q<1l) 
11  mn 


-  2Rd,,A  mnt  sintmt^costptOcostntWsi nta*1!) 

l  o  mn 

-  R2d,.A  n2t2cos<nrt*)co*<ptt)*in<nt|l)*in<qt'S) 

6g  mn 

-  dj^Brwim2<2sin<mtt)co»<pt4)co»(nt1))sin(qt11> 

-  Rtdlo+d,,)B  ran^costmt^costpl^sintntDsintq*1)) 

l  /  oo  mn 

-  R2d_ .  B  n2t2s  i  n  <m<%)  c ost  pt%)  cost  n*1))  s  i  n  t  qt1)) 

110  mn 

-  ks2»g5Amncostnrt(i)costpt%>sintnt1))si  ntqt1!) 

-  k sa^C^nrtc ost mil)  cost  pt%)  *  i  n  ( nit1))  s  i  n  t  qt1)) 


♦  («  /12)Afancot(mlC%)co»(pUl)*in(nlt'rl)sin(qtT))  d%  dl 


IRAmnm(cos<nt1))*in(qf1l>  *■  Bronm*costr»1(1))5intq*1l> ) 


-  dj^[RAmnntco*<mt) cost p<) cost ntDs-intqt1!) 


4  Bmnmtco*<m1t)c04<plC>co*<n1t1>>5in<(Tt1'> J  dl  ■  8 


v>-:v 
%  *• 1 


a 


From  Eq.  (76) 


Ziff  -d,  .A  n2t2cos<mtl>*in(pt|)sin<nt1l)cos<qt1l) 

r-ln-lJeJg  16  m 


-  R<d,„+d..)A  mnlt2si  n<mt%)»in<p<4)cos<ntl)cos<qlC*l) 

l  z  66  ion 

-  R2d_.A  n2t2co»(mt%)*in(pt4)sin<nt1))co*<qt1») 

Z6  mn 

-  d..B  ro2t2*in(rrt%)»in<pt%)co»<n<1l)co»<q<1i) 

66  mn 

-  2Rd_  . B  mnt2co*(rM(ll)sin<plC%)sin<nl*l)cos<q‘<‘n) 

do  mn 

-  R2d0_B  n2<2»in<mlC%)sin<pt*)cos<n<il>cos<qt*> 

ZZ  mn 


-  k*  *44Bmnfcin<nrtC*l)*in<pt%)co»<nt'H)co*<qtl> 

-  kRao^C^nts  i  n  <mt%> » i  n ( pK%> cos<  ntl)  co*(  qtrt) 

+  (w2/12)B  *in(ml%)sin<pH)co»(n<Tl)co*<qt1J)  dl  d* 
mn 


f 

•'a 


*  I  d_.IRA  ntcos(ml%)tin(pl4)  *  B  mlcos<nl%)si n<pt%) ] 
I  g  do  mn  mn 


-  d_,[RA  nfco»<nt)cos<qf)cos<mt%)sin(ptt) 
Z6  mn 


♦  B(tnmfcoa(nt)cos<ql)co*<ml%}sin<p<|)3  d*  ■  8 


<73) 
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When  m  *  p  (or  n  *  q) 


cos(mtx)  cos(ptx)  dx 


sin(nrttx)  sin(ptx)  dx 


sin(mlx)  cos(ptx)  dx 


sin(m-p)tx 

2t(m-p) 


sin(m-o)tx 

2t(m-p> 


cot(m-p)l(x 

2t<m-p> 


sin(m-»p)tx  1* 
2t(m+p>  |0 


sj_n(m±Eii2. 1 1 

2t(ra+p)  |0 


cos<mtp)^x  1* 
2t<tn+p)  I 


0 


0 


cos(ra-p)*  -  1  _  cos(m+p)t  -  1 
21(m-p)  2t(m+p> 


<-cos<ra-p)<  +  i)<m+p)  +  (-cos(n>+p)t  *  l)(n>-p) 
2tC(m-p)<m+p)  21(ra-p)(rn+p) 


<-cos(n>-p)t  +  lXm+p)  +  (-cos<ro+p)t  »  l)(nrp> 

2  2 
2<m  -  p*> 


■  •  ■for  (m+p)  an  •  v*n  integer 


(88) 


and 


2  2 

■  2m/*(m  -  p  )  for  (m+p)  an  odd  integer 


(81) 


Similarly  for 


s: 


tin(pKx)  coa(mKx)  dx  ■  0  for  (p+m)  an  avan  intagar 


2  2 

»  2p/K(p  -  m  )  for  (p+m)  an  odd  intagar 


(82) 


(83) 


Wa  ahall  now  introduca  tha  notation 


cotimlCx)  coa(ptCx)  dx  •  Oi  or  I) 


(84) 


<  78) 


<7 9) 


where  tha  valut  on  tha  laft  aida  of  tha  paranthalical  expression  it  tha  va’ue  of  tha 
integral  whan  m*p  and  tha  value  on  tha  right  side  of  tha  expression  is  whan  m%p. 


Thu*,  upon  integrating,  simplifying,  and  taking  only  a  finite  number  of  terms  we 
have  (note:  even  and  odd  tests  are  for  the  *um  of  the  indices): 


From  Eq.  (72) 


M  N 

2  X  C t204  or  0)04  or  0)<-m2dn  -  n2R2d46  "  ks2a55/«2  +  S2/12t2) 

nif*!  n=l  ' 


2  2  2  2 

<0  or  0  even,2n  odd)<0  or  0  even,2q  odd) <2Rmndj^)/< (m  -p  Xq  -n  )) 


2  2 

+  (0  or  0  even,2q  odd)  <nRd^<  l-cos(mt)  cos(p<) )  )/(q  -n  )]  Amn 


+  It  <X  or  0)04  or  0)<-mnR<d12  +  d^,)) 


-  <0or0  even, 2m  odd)(0or0  even,2q  odd) <m2dj^+n2R2d9/.)/( <m2-p2) (q2-n2) ) 


26 


2  .2, 


♦  (0  or  0  even,2q  oddXmdj^d-  cos(mt) cos(pt) ) )/(q  -n  )]  B( 


mn 


-  [104  or  0)04  or  0)(mksa„)l  C 

jd  mn 


XX  0 


(85) 


From  Eq.  (73) 


M  N 


X  X  tt‘04  or  0)04  or  0X-ranR(dlo  +  d,,)) 
,  .  12  66 
rw=l  n=l 


-  (0or0  even,2p  odd)(0or0  even,2n  odd)(m2d1(j+n2R2d2^>/<<p2-m2)(n2-q2)) 


2 _ 2. 


MB  or  0  even,2p  odd) (nRd_ ,d-cos<n*)cos(qlC) ) )/<p  -m  )]  A 

<LO  l 


mn 


♦  C1204  or  0) 04  or  0Xto2d6|S  -  n2R2d22  -  (ks2a44>/t2  +  <32/12t2) 


2  2  2  2 

<0  or  0  *ven,2p  odd)(0  or  0  even,2n  odd)<2Rmnd2^)/<(p  -m  Xn  -q  )) 


♦  <0  or  0  *ven,2p  oddXrad»,< l-cos(nt)cos(q<) ) )/(p2-m2)  1  ii 

mn 


£104  or  0)04  or  0) <nRks*.. ) J  C 


(86) 


And  from  Bq.  (74) 


M  N 

Z  Z  -  <04  or  0)04  or  0)<mksacl->  A 
.  .  55  ran 

iw*l  n=l 

-  <04  or  0)04  or  0)<nRksa„il>  B 

44  mn 

+  [l2  04  or  0)04  or  8)(-n2ka„-n2R2ka.. 

55  44 

♦  k  2\ j  m2  <  R2/ s2 ) + k  j  X  j  n  2  <  R4/ »2+x2  <  R4/!2  s2 ) ) 

♦  <0  or  0  even,2p  odd)<0  or  0  even,2q  odd) 

♦  <2kgmnR\j)/<52<p2-m2Xq2-n2) )  ]  C  «  0  (87) 

Equations  (85)»  186),  and  (87)  are  now  reedy  to  ba  programmad.  Th»  process 
to  generate  the  Galerkin  equations  is  to  pick  an  M  and  N  (they  must  be  equal),  then 
cycle  through  p  and  q  (which  are  equal  to  M  and  N).  Thus,  M  and  N  determine  the 
number  of  terms  in  each  Galerkin  equation  and  p  and  q  determine  the  number  of 
Galerkin  equations. 


For  a  plate  clamped  on  all  side*(  on*  Know*  from  plate  theory  that 


In  order  to  substitute  into  Eqs.  (56),  (59),  and  (63)  one  needs  the  following 


derivatives  (note:  all  ES  are  over  m=i  to  •  and  n=l  to  •>. 


2 

2 

A 

mn 
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= 
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= 
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2 
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mn 
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s 

2 

2 

B 
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2 

2 

-B 
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2 

2 

C 

ran 

w,« 

x 

2 

2 

-C 

ran 

w,^ 

x 

2 

2 

C 

ran 

w,i 

X 

2 

2 

C 

ran 

w  __ 

X 

2 

2 

-C 

mt  cos(mt$>  sinfnt1)) 

m2t2  sin(rati)  sinCnt*)) 

rant2  cos(ntt)  cos(nti)) 

nl  sin(nrt$>  cosfnl1!) 

n2t2  sin(rat%)  sinint1!) 

rat  cos(m t$>  sin<nt“9) 

ra2t2  sin(nrt|)  sin(nfn) 
2 

rant  cosintfc)  cosint1)) 

nt  sin(mt^)  cosCnt1!) 

nV  sin(ret%)  sin(nfn) 

rat  cos<mt^>  sinCnt1)) 

ra2t2  sln(ntt)  sin(nttl) 
2 

rant  cos(ratt)  cos(ntTl) 
nt  sin(mt^)  cosfnt1!) 
n2t2  sin(ratt)  sinint1)) 


(92) 


Now,  substituting  into  Bq.  (36)  w*  havt 


I  iff  m2l2sin<mtq)sin<n<n)  ♦  2Rd.  ,A  mn*2cos(n»t%>cos<nt't> 

n*l*W8  11  mn  16  mn 


-  R2d.,A  n2t2sin<«nt^)sin<n<1l)  -  d,  ,B  r»2t2sin(mH)sin<nt1)) 

66  mn  16  mn 

+  R(d. ^+d  )B  rant2cos<mt*)cos<nl1>  -  R^.^B  n2<2sin<mt*)sin<n-n> 
i  2  oo  ran  26  mn 

2 

-  ks  sin(mt^)sin<nt1t)  -  ks*cl-C  ra<cos(ra<%)sin(nt'n> 

do  mn  55  nn 


+  <u  /12)A  sin<mti)sin<nt^))  (sin(pU)sin(qtn)}  d*  dl 

nn 


+  I  <  d.  .A  mtcos<8)sin<nf^)  +  Rd,-B  nfsin(8)cos(nf11)) 
Jq  1 1  mn  1 2  mn 


♦  dJ(i[RAmnntsi ntBicostnl^)  ♦  BmnrsKcos(0)sin<niC,l)])  (sinCBJsintqt:11))  > 


{  d,.A  mtcos<mt)sin<n<1t)  4  Rd.^B  ntsintnnDcosCntR) 
1 1  mn  1 2  ran 


+  d-.IfW  ntsin<mt)cos<ntt)  ♦  B  nrtcostmDsintnl11))  J} 
1 6  mn  mn 


*  <sin(pt)sin(qt'll) )  d1! 


r< 


.A  nKcos<ratcq)si n<8)  ♦  Rd-.B  nts i n(ntf%)cos<0) 
6  ran  26  ran 


4  d,.TRA  ntsin(mf^)cos<0)  ♦  B  ra*co5<n<%)sin<8>3>  (sin<p<*)sin<0)} 
66  mn  mn 

-  (  d.  .A  mtcos<nrt%)sin<nl)  4  Rd» .B,ntsin<nflC%)cos<nt) 

16  mn  26  ran 

4  d,,[RA  ntsin<rat%)cos<nt)  4  B  mtcos<mt%)sin<nl(>]} 

66  mn  mn 


*  (si n<ptq>s j n(qt> )  d%  ■  0 


(93) 


Substituting  into  Sq.  (59)  we  have 


S  I  (Tt-d 

n^ln*lwW0 


,.A  m2t2sin(ra<%)sin<nt‘i))*R<d 
16  nn 


2 

.  ^d.,)A  mnl  cos<ml%>cos(nlil) 
l c  66  nn 


-  R2d,  A  n2l2sin<ml$)sin<nt4)  -  d..B  ra2l2sin<m«>sin<nli» 

Li*  run  66  mn 

♦  2Rd0,B  mnt2cos<ml$)cos<nH)  -  R2d0_B_  n2l2sin<mlll)sin<nl‘1)) 

26  nn  22  nn 

2 

-  ks  *..B  tin<ral%)sin<nll)  -  kRsa-.C  nisi n<m!4) cos (nlH) 

46  ran  44  ran 

♦  <52/12)B  sin<mH>sin<nl*l)}  (sin<pH>sin<qli))}  d*  dH 


mtcos(8)si  n<nH)  ♦  Rd 


26^mn 


nlsin(0)cos(nl1)) 


♦  d,.[RA  nlsin<0)cos(nl1))  +  B  mlcos<0)sin(nti))]}  <sin(0)sin(ql1)>} 

66  ron  mn 

-  (  d.  .A  nrlcos<ml)si n<nl1))  +  Rd~,B  nlsinCmDcosCnl1!) 

16  mn  26  ran 

+  d  [RA  nisi  n<ral)  cost  nil)  ♦  B  mlcostmDsintnl11!)]} 

66  ran  mn 

*  (si n<pl)si n<ql1)) )  d4 


d1oAimiiralcos(ral%)sin<0)  +  Rd_„B.nlsi  n(mll>cos<0) 
i  i  nn  22  nn 


♦  d_.[RA^  nlsin<raH)cos<8)  ♦  B  „nr!co*<nH)sin<0)]}  <sin<pl*)sin<0>) 

l  o  mn  mn 

-  t  d,_A  mlcos<mH)sin<nl)  ♦  Rd„»B  nlsin(ml|)cos<nl) 

12  mn  22  ran 

♦  d^.tRA  nlsin<n»H)cos(nl)  ♦  B_  mlcos(nl%)si n(nl) ] } 

26  nn  mn 


*  (sin(pli)sin(ql) }  d(  >  8 


<94) 


And  finely  *ubstituting  into  Sq.  (63)  we  have 


Z  Z  |* J*  £  k5a5^mnmtcos<mt%)ftin<n<T))  -  kag5Crerim2K2sin(ral^)sin(n<1)) 
rtf*  1 

-  kRsa.^B  ntsin(mt%>cos<ntH)  -  kR2atf4C  n2t2sin(n<^)sin<nl1)) 

44  tin  mn 

♦  k^j  (R2/*2)Cmnm2t2sin<rol4)sin<n<1)) 

♦  2k„(R3/s2)X,C  mn<2cos(nrKOco*(nt‘t) 

3  l  mn 

+  k.(R4/*2)x,C  n^sintnrtDainint1)) 

1  l  mn 

♦  <R4/s2)x.C  »in(m<t)»in<nt1))}  {sin<p<S)*in<qt1))>  d%  d1) 

2  mn 

♦  r1{ksa5^mnsin<0)sin<n<i»  ♦  kaggC^mlcosCSJai n(nti)) )  Uin(0)sin(qtT))3 
8 

-  {ksag5Aransin<mt)sin(ntH)  +  kaggC^ertCcoadiOsi n(nt^l) ) 

♦  <sin<pt)*in<q<‘4)}  dl 


r 


(k*a.  .B^sintmUJsiniB) 
44  etn 


Rka. .C  nt*in(nrlC$.>co*<0)}C*in(p<%>*in<0)5 
44  mn 


-  (kaa .  .B  sin(mt4)*in<nt)  +  Rka.-C  ntsi n<*rit%>cos<n<> > 

44  mn  44  ran 

*  (sin<pt%)sin(qt)>  dt  ■  0 

The  indicated  multiplication*  are  performed  after  noting  that  sin(0)  B 
*in(p«)  ■  *in(qi()  ■  0  (p  and  q  are  integer*)  and  cos(0)  *  1. 


From  Sq.  (93) 


•  -  /»1  -1 

nil  n-lJeJe  ~dllAmnt  *in<at*>tin<Pt»sin 


<nt1))sin<qt'i)) 


+  2Rdi4Amn™nt  c os (m*4 ) » i n < ptt ) c o*< nf1)) s i n < qfq) 

2  2  2 

R  d66Amnn  *  4in<rat')*in<P<%>*tn<n<^)»in(qt^> 

2  2 

dl6Bran™  *  sln<™*4>»in<ptq>»jn(nti»>»jn(qH)) 

2 

R<d12+d66)Bmnrant  co*<m<*>*i rXplOcostntDsiiUqtH) 
2  2  2 

'  R  d26Bmnn  1  tin<mt^>B'n<Pt%>*in<n*i))*in(q«) 


- 


*55^1*0  * ' n  <n**>  *' « <  P<4>  si  n  ( n<*))  %  i  n  ( qlH) 


“  K**55C»nm<c os<rrt% ) s i n  <  pt*>  * i n ( n«H) %  I n ( qt n> 


,a2 


♦  <U  /12)Amn*in<mt^)*in<pt%)*in(ntH)*in(qtH)  d%  d4  »  0 


(96) 
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From  Eq.  (95) 


I 

n* 


i  •'e  •'e 


ksa__A  m*cos(mt^)sin<p*$)  si  n(nO)si  n(qt'n) 
D5  run 


-  k a^C  m2t^sin(ml?)sin(pl(%)*in(nt',))sin<q<^) 

dd  mn 

♦  kR5a44Bmnn<5in<mt^)sin<p<4)co«-<nt1,l)sin(qt‘n) 

-  kR2a,.C  n2l2sin(rrti)*in<pt%)sin(nl1))sin(q<1!)) 

44  mn 

+k2x1<R2/*2)Cmnm2t2sin<ml%)sin<pt^)sin<n<il)sin(qiCil) 

+  2k0(R3/s2)x.C  mnl2co*(m<4)5in<pl%)cOB(nl(‘n)sin<ql‘9) 

J  1  mn 

+  k,  (R4/s2)X.C  n2t2*in<mt%)*in(p<%>sin(nlCil)si  n(q<i|) 
i  l  mn 

+  (R4/»2)x„C  sin<nrt(%)*in<pt%)*in(ntt1)*in<q<1l)  *  0  (98) 

i  ran 

We  can  now  integrate  Bqs.  (96),  (97),  and  (98)  using  the  same  integrals  from  the 
previous  section.  Thus,  upon  integrating  ,  simplifying,  and  taking  only  a  finite 


number  of  terms  we  have 


From  Eq.  (96) 


M  N 

Z  Z  It204  or  0)04  or  •><-«2d1]  -  n2R2d^  -  ks2»55/<2  + 

m“l  n=l 

-  <0  or  8  *y»n ,2p  odd)(0  or  0  *v*n,2q  odd) 

*  (2Rmnd,  .)/( <p2-ra2Xq2-n2) )  ]  k 

16  mn 

4  It2Oi  or  0)04  or  0)<-t»2d16-n2R2d26) 

4  (0or0  even,2p  odd)<0or9  tv*n,2q  odd) 

#  (mnRtd.^+d,  .))/<< p2-m2)<q2-n2))l  B  ^ 

1 L.  66  IBn 

-  ll<8  or  0  even,2p  odd)  04  or  0)<mk*i^g)/(p2-ra2) )  Cffln  *  0 


From  Eq.  (97) 

M  N  2  2  2  2 

Z  Z  [<*04  or  0)04  or  0)(Vd,  ,-n*R *d„.) 

.  .  16  /d 

i*»l  n*l 


-  <0or0  evcn,2p  odd)(0or0  evtn,2q  odd) 


*  (ranR(dto4d.  .))/<<p2-tn2Xq2-n2))] 

1 4  66  ran 

4  t<204  or  0XH  or  0X-ia2d(S6  -  n2R2d22  -  <ks2*,44>/<2  4  52/l2t2) 


4  (0  or  0  tv»n,2p  odd)<0  or  0  ev«n,2q  odd) 


#  <2ft*nd_.)/<<p2-m2Xq2-n2)>3  B 
26  mn 


2  .2, 


-  [<<*  or  0X0  or  0  *w*n,2q  odd)<nRks*44)/<q  -n  >J  C 


ran 


-2  2 
(i)/12«) 


<??) 


(100) 


And  from  Eq.  (98) 


2  2 

2  S  -  t<0  or  0  euen,2p  odd)04  or  0)<mksa__)/<p  -m  >  A 
m=l  n-1  55 

♦  <04  or  0X0  or  0  *w*n,2q  oddXnRk*a..)/(q2-n  )  B 

44  ^  ran 

♦  t<2  04  or  8)04  or  8)<-fli2l(a„-n2R2ki„ 

55  44 

♦  k2Xjm2<  R2/s2) +k  Xj n2( R4/s2+X2 ( R4/l2s2 ) ) 

+  <0  or  0  ev»n,2p  oddX0  or  0  even,2q  odd) 

♦  <2k^nnR3x1)/<s2(p2-ra2Xq2-n2))]  Cmn  =  0  (101) 

Equations  (99),  (100),  and  (101)  art  now  ready  to  be  programmed  to  generate 
the  GalerKin  equations  as  before. 


52 


O  '  'l*.'  ^  "  . 


--  •-  -- 


For  a  plat*  clamped  on  two  opposite  tides  and  simply-supported  on  two 


opposite  sides>  one  has  the  following  boundary  conditions. 


at  x  ■  0,a 


<i«2> 


and 


w  ■  f  ■  »  *  0 

x  y 


at  y  *  0,b 


w  »  ?  ■  0 

x 


(103) 


+  D_,  f  +D0,  <T  +  <r  )  *  0 
y  12  x»x  22  y*y  26  y,x  x,y 


Therefore,  w*  shall  choose  for  our  admissible  functions 


f  -  Z  Z  A  sin<mtx/a)  *in(nty/b> 
x  <  _  i  ran 

iw*l  n-1 


*  -  I  Z  EL  sin<rnlx/a)  cos<ni(y/b) 

y  «w«l  n-1  mn 


W  -  Z  Z  C  sin<mtx/a)  sin(nlCy/b) 
m*l  n-1  mn 


or  in  terms  of  normalized  coordinates 


L  —  Z  Z  A  sin(raK^)  sin(nl^) 
*  w-1  n-1  mn 


Z  Z  El  tin(mft)  coslnl1!) 
ra-1  n-1  mn 


W  -  Z  Z  C  sin(ra<%>  sin(n<^l) 
wl  n-1  mn 


(104) 


(105) 


To  substitute  into  Bqs.  (5 6),  (59),  and  (63)  we  need  the  following  derivatives 


(note:  ell  ££  *re  over  m*i  to  «*  and  n«i  to  *>. 


3 

£ 

£ 

A 

an 

a 

£ 

£ 

-A 

an 

3 

£ 

£ 

A 

an 

Y* 

m 

£ 

£ 

A 

an 

«* 

Yvl 

m 

£ 

£ 

-A 

an 

Y* 

m 

£ 

£ 

Bran 

Y*% 

m 

£ 

£ 

-Bmn 

Yv) 

3 

£ 

£ 

Bmn 

Yl 

3 

£ 

£ 

”Bmn 

Y*1* 

m 

£ 

£ 

-B 

an 

w,% 

8 

£ 

£ 

C 

an 

W,M 

m 

£ 

£ 

c 

oe 

1 

w,%1 

m 

£ 

£ 

Swi 

w,' 

s 

£ 

£ 

Can 

5  .. 

. 

£ 

£ 

-C _ 

ret  cos(ittH)  sin<nt1)) 
re2!2  sin(mt%)  sin(nt^) 
mnt2  cos(nrt%)  cosfnt*)) 
nt  sin<mt%>  cos(ntt) 
nVsi  n(at4>  sin(nti)) 
at  coiim&V  cos(ntt) 
mV  sin(at%)  cos(nt*l) 
mnt2  cos(ati)  sin(nt*l) 
nt  si  n<nrt%)  sin(ntt) 
nV  sin(ati)  costnt*)) 
at  cos(mt%>  sin<ntt) 
mV  sin(att)  sin(nt1l) 
mnt2  cos(at%)  ccsfnt1)) 
nt  sin<mt|>  cosfnt1)) 
nVsi  n(atl)  sin(nt^) 


<  106) 


Now  substituting  into  Eq.  (36)  w •  havs 


2 

ra«*l 


2  Pf  <-d,4A  a2t2si  n<rat%)sin<n<i))  ♦  2Rd,  .A  mnt2cos(m<%>cos<nt't> 

n-lJeJe  11  ron  16  mn 


-  R2d,.A  n^t^sin<s»H)*in(nt1))  -  d.  .B  a2t2sin<mt$)cos<n<1,l> 

66  mn  1  o  mn 

-  R(d.0*d,.)B  mnt2cos<mt$)sin<nt*l)  -  R^d^.B  n2<2sin<mt%)cos(nt‘,l) 

12  66  an  26  mn 

-  ks2a__A  sin(nrit%)sln<nl1))  -  ksa__C  m<cos(rat*>si n<nl1)) 

mn  55  mn 

♦  (G2/12)A  sin<mf%)sin(nl1))}  {sin(ptOsin<q1tn)}  d*  dH 


d,  ,A 
11  mn 


mlcos<8)sin<nl‘l)  -  Rd,_B  ntsin(8)sin(nt1l) 
i  l  i  mn 


♦  d.  .[RA  ntsin^costn*1))  ♦  B  nrlcostaJco^n*1)}])  {sin(8>sin<qt^) ) 

1 6  mn  mn 

-  (  d,  ,A  rotco»(nrt()sin<nt1,l)  -  Rd.-B  ntsin(niK)sin(ni(’)) 

1 1  mn  1 2  mn 

♦  d,.[RA  n‘tsin<mt)cos<nlC1)  *  B  mtco*<mt)cos<nl’l>3} 

16  mn  mn 

♦  {sin<pt)sin(ql1))}  dl 


{  d,  .A  afcos(af%)sin<8)  -  Rd„,B  nt*in(nrK%)sin<8) 
16  an  26  an 


+  d,.tRA  nfsin<mft)cos<8)  ♦  B  m<cos(nt%)cos<8) ] )  <sin(pK^)sin<0)> 
66  mn  mn 

-  <  d, .A  mtcos<mti)sin(n<)  -  Rd0-B  ntsin<mft>sin<nt) 

1 6  mn  26  mn 

♦  d..tRA_  nisi n(mt%)cos(nl)  +  B^  mtcos<mt%)cos(nlO 3 ) 

66  mn  mn 

*  (*in<pt%)sin(qt)3  d*  ■  8  (18?) 


55 


Substituting  into  Bq.  (59)  we  have 


•  •  22  2 
Z  Z  \  \  {-d.  .A  m  i£sin(mt^)si  n<n<10  +  R< d.  _+d  .  ,)A  rrrnt  cos<mt^)cos(ntT)) 
linn  l  2  44  mn 

m=l n=lW8W0 


-  R^d0.A  n^t^sin<srtt%)sin(nttl)  -  d.,B  ra^t^si n<mt%)cos<n<1l) 
24  ran  44  ran 


2Rdoz8  mn<2cos(nrtC%)sln<nt1l)  -  R^d_„B  n^t^sin<mt§)cos<rvr,l) 
26  mn  22  mn 


ks  a..B  sin<nrC%)cos<nt1l)  -  kRsa-.C  n-ts^mt^cosCn*1)) 
44  mn  44  mn 


+  (5V12)Bi>nsin(nrt|)cos(n<i)))  <sin(p<*>cos<q<i))}  d%  d^l 


J/  ^ 


mtcos(0)si n<nfl)  -  Rd~,B  nfsin(0)sin(nt'n) 
mn  26  mn 


♦  d .  IRA  ntsi  n(0)cos(nt1))  ♦  B  nKcos(0>cos(nf  ^l)  ]  }  {si  n(0)co*<qtH)  } 

66  mn  mn 

-  (  d.  .A  mfcos<mrC)sin<nf1))  -  Rd^.B  mtsin<m<)sin(nl(i)) 

16  mn  2o  mn 

♦  d,,[RA  ntsin<mrK)cos<nt,n)  +  B  mlcos<mt)cos(ntTi:)  ] ) 

66  mn  mn 

♦  (si n<pt)cos<qt1)) }  d1! 


mtcos(nrit$)si n<0) 


-  Rd-^B  n<sin(m<4)sin<0) 
22  mn 


♦  d_.[RA  n<sin(nrt(%)cos<0)  ♦  B  mfcos(mlC%)cos(0>]>  <sin(p*t)cos<0) ) 
26  mn  mn 

-  (  d. 0A_vi%mtcos<mt^)sin(nt)  -  Rd-_B  ntsin(rat*)sin<n<> 

12  mn  22  mn 

+  d,.[RA  ntsin<mt%>co*<nt>  +  B  mtcos<tnt4)cos<ni(>  ]  } 

26  mn  mn 


•  (si n(pl%)cos(qf ) }  d%  ■  0 


(108) 


And  finally  substituting  into  Bq.  (63)  we  have 


-  (Ksa. .B^sin(mH)cos<n<)  ♦  RKa^-C  ntsi n(rat%)cos(nf) } 

44  ran  44  ran 

*  {sin(pt*)sin(qiO)  d%  *  0  <10?) 

Again,  performing  the  indicated  multiplications  and  noting  that  sin(0)  ■ 
sin(pt)  «  sin(q*)  ■  0  (p  and  q  are  integers)  and  cos<0)  ■  1  we  have 


From  Eq.  (108) 


Z  E  f  (*  -d,  A  ra2l2sin(mt$>sin<plC%)sin<nlC1ll)cos<qt‘,l> 

vl  n-l  16rmn 


+  R(d.0+d..)A  mn<  cos<mt^)5in<pt%)co»<nl1))cos<qK1)) 
l  x  66  ran 

-  R2d„.A  n2*2air»<nrt£%)ain<p<l)ain<n<1l>coa<q't‘,l> 

26  inn 

-  d.,B  m^t^*in<mt%)*in<pl^)co5(nt1l)co*(ql1)) 

66  ran 


-  2Rd„,B  mnt  co»<iMt%)*in<plC%)sin(nl(1l)co*(ql1,l) 

26  inn 

-  R2d00B  n2t2sin(m<%)ain<pl(*}cos<nf*))coa<qt1l) 

22  inn 


-  ks  a44Bransin<ml%)si n<pt%)co»<nl1[))co»<q'C1l) 

-  kRaa^C  nta  i  n  ( mt%>  a  i  n  (  pt%) c  oa(  nH>  coa<  qt1)) 


♦  («  /12)B  sin<mt%)sin(pt%)cos(nt1))cos(qfi))  d$  dl 

ran 


Jd,,[RA  ntnin(mH)sin(p*%)  ♦  B  mKcos(rat^)si n(pl%> ] 
0  x6  mn  mn 


d  - .  [  RA  ntc  oa<  nt ) c  os  <  q* ) a i n  (inf % )  a  i  n  <  pt%) 
26  ran 


♦  Bnnrrtcoa(n< )coa(qt)coa<mtl)sln(p<$) ]  d%  ■  0 


(111) 
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From  Bq.  (169) 


2  2 
w»l  n*l 


nrttc  os(rat|) » i  n  <  ptl>  *  i  n  ( nt^l)  s  i  n  ( qtH) 


-  ka«C  ra2<2sin(mt%)sin(ptO*in(nH)sin(qt‘*)> 

jj  mn 

-  kRsa. .B  ntsin(ml%)sin(pt$)sin(nt1l>*in(q1Ci)) 

44  mn 

-  kR2a. .C  n2l2s i n (ra<%) s i n < pt|) s i n ( nt^l) s i n ( qt"1)) 

44  mn 

♦  <  R2/s2  )  Cmnm2<2*  i  n  (mt* )  s  i  n  ( p«* )  3  i  n  <  nil)  *  i  n  ( q<D) 

+  2k,,(R3/s2)x,  C  mnt2co»(mt4)sin<pl%>co*(nl1))*in<qt1)) 

J  1  mn 

♦  k,<R4/s2)X,C  n2t2sin(mt%)sin(ptl)sin<nin)sin(q'C1l> 

1  1  mn 

♦  (R4/s2)x_C  sin(mrt|)sin(p*%)*in(n<1))sin<qt1l>  cHUrt  »  8  (112 

c  mn 

We  can  now  integrate  Bq*.  (1 10)*  (1 1  !)■  and  (112)  u»ing  th*  tarn*  integrals  from  the 
previous  sections.  Thus  upon  integrating!  simplifying!  and  taking  only  a  finite 


number  of  terms  we  have 


From  Sq.  (114) 


H  N  o 

Z  Z  ll2<*  or  ex*  or  0)<-m2dH  -  n2R2d^  -  ks2a55/*2  +  W2/12*2) 

rr«l  n=l 


Ml  or  0  «v*n,2p  odd)(0  or  0  evtn,2q  odd) 


#  (2Rmnd,  .)/<  (p2-n2Xq2-n2) )  1  A 

16  mn 


2 _ 2S 


-  l«0  or  8  »v*n,2p  oddX*  or  eXmnRCd^  4  d^))/(p  "m 


-  t(*  or  0X4or§  tv*n,2q  odd)<m2dj^+n2R2d,A>/<q2-n2) 3  B 


26 


mn 


2  _2. 


-  [(I  or  1  *v*n,2p  oddX*  or  0Xmk»a^)/(p  in  >3  C 


mn 


(113) 


From  Bq.  (Ill) 

h  N  -  2 

Z  Z  I«0  or  0  *v*n,2p  oddX*  or  0XmnR(d12  ■»  d(S6>)/<p  -m  ) 
ro*l  n=l  * 

-  t(*  or  0X0or0  #v»n,2n  oddXi»2dj£*n2R2d2£>/<n2-q2> 

♦  t<*  or  0XnRd_.<l-co*(nt)co»<q<)))3  A 

26  mn 

4  [<2<*  or  0 ) <*  or  0)<-m2d, ,  -  n2R2d,_  -  (ks2a..)/t2  ♦  Q2/12t2) 

66  22  44 

-  (0  or  0  »v»n,2p  oddX0  or  0  *v«n,2n  oddX2Rmnd2£)/((p2-m2Xn2-q2) ) 

♦  (0  or  0  #v*n,2p  oddXmd, .< l-co*(nt)co*<q<) ) )/<p2-m2)  1  B 

£u  mn 

-  [«*  or  0X*  or  0)<nRksjL,.)]  C  -  0  (114) 


And  from  Sq.  (112) 


M  N  _  „ 

Z  Z  (8  or  0  eM«n,2p  odd)  04  or  0)(mksa_..)/(p£-flr )  A 
4  4  n 

iw«l  n»l 


-  <04  or  8) (H  or  0)(nRksa..) 

44  mn 


♦  It2  04  or  8)04  or  0)(-ra2kaee-n2R2ka, 


+  k2x1m2<R2/*2)+k1x1n2(R4/*2+x2(R4/t2»2)) 


♦  (0  or  0  even,2p  odd) (8  or  0  even,2q  odd) 


*  <2k^nnR3X1)/<*2(p2-n2Xq2-n2))J  Cn 


(115) 


Bquations  (U3)i  (114),  and  (115)  art  now  ready  to  be  programmed  to 
generate  the  Galerkin  equations  as  before. 


.'V  N* 
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This  chapter  will  describe  the  computer  program*  used  to  solve  the 
eigenvalue  problem  formulated  in  the  last  chapter.  It  will  also  discuss  the  two 
types  of  laminated  plates  used  in  the  programs  and  the  subsequent  analysis 
performed  with  those  plates.  We  will  begin  by  describing  the  computer  algorithms. 

Computer  Algorithms 

Four  computer  programs  were  written  to  set  up  and  solve  for  the 
approximate  natural  frequencies  and  mode  shapes  for  the  three  boundary 
conditions  considered.  The  first  program  determines  the  nondimensional  bending 
and  extensional  stiffness  elements  for  a  symmetric  laminate.  The  second  program 
formulates  the  eigenvalue  problem  CA3xaXCB3x  where  CA3  is  the  stiffness  matrix 
and  CE3  is  the  combined  mass  and  inertia  matrix.  The  third  program  solves  the 
eigenvalue  problem.  The  last  program  generates  the  data  for  the  mode  shape 
plots.  A  description  of  each  program  follows. 

Program  One  determines  the  nondimensional  stiffness  elements  and  the 
density  of  the  plate.  It  is  divided  into  three  basic  sections  and  allows  for  a 
laminate  to  be  built  up  from  different  materials  and  orientations.  The  three  basic 
sections  are  1)  the  input*  2)  the  computation  of  the  dimensional  bending  and 
extensional  stiffnesses*  and  3)  the  computation  of  the  nondimensional  stiffnesses. 
The  input  section  gathers  the  following  data: 


1)  Orientation  angle*  © 


Poisson's  ratio*  Vgj,  it  calculated  from  by  the  tquation  x  (Eg/ E^). 

The  shear  moduli  terms*  Ggg  and  Gigt  are  computed  from  G^.  GJ3  it  set  equal  to 
Gj2  and  Ggg  is  set  equal  to  88X  of  Gjg. 

The  second  section  of  this  program  is  the  longest.  It  calculates  the 
dimensional  bending  and  extensional  stiffnesses  for  a  lamina  and  then  loops  back 
to  the  beginning  for  more  input  data.  If  the  material  properties  are  the  same*  it 
will  only  ask  for  the  new  orientation  angle.  It  will  continue  to  sum  these 
stiffnesses  together  until  the  laminate  lay-up  is  completed.  The  stiffness 
elements  are  determined  by  computing  the  reduced  stiffness  terms  (Q^'s)  from  Bq. 
(14)  and  then  computing  the  transformed  reduced  stiffness  terms  (Q^'s)  using  Eq. 
(13).  The  extensional  and  bending  stiffnesses  are  then  determined  using  Eqs.  (28) 
and  (3t).  After  the  laminate  buildup  is  complete*  the  program  moves  on  to  the 
third  section. 


Jk 

•*  . . 


The  third  section  takes  the  dimensional  stiffnesses  and  computes  the 


nondimensional  stiffnesses  using  Bq.  (53)  from  the  previous  chapter.  These 
stiffnesses  are  then  printed.  A  listing  of  this  program  may  be  found  in  Appendix 


£  > 

'  I 


The  second  program  computes  the  mass  and  stiffness  matrices  for  the 
eigenvalue  problem.  Like  the  first,  it  too  is  divided  into  three  sections.  The  first 
section  contains  the  plate  data  generated  from  Program  One.  The  second  section 
builds  the  stiffness  and  mass  matrices  and  the  third  section  outputs  these 


matrices  to  a  file. 
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Section  one  is  the  input  section.  It  contains  all  o-f  the  required  material 
properties!  geometrical  plate  data*  and  integers  which  determine  the  number  o f 
terms  in  GalerKin's  equations  and  the  number  of  Oalerkin  equations. 

Section  two  generates  the  stiffness  and  mass  matrices.  It  takes  the 
algebraic  equations  (Bqs.  85-87  for  the  first  boundary  condition*  Bqs.  99-18 i  for 
the  second  boundary  condition!  and  Bqs.  113-115  for  the  third  boundary  condition)* 
determines  the  value  of  the  integrated  sine  and  cosine  terms*  and  then  builds  the 
family  of  Oalerkin  equations.  The  stiffness  and  mass  (or  inertia)  terms  are 
separated  and  the  matrices  are  constructed.  It  should  be  noted  here  that  the 
buckling  terms  (those  multiplied  by  kj,  k^i  and  kg)  were  not  used  in  this  analysis 
and  will  not  appear  in  the  computer  listing. 

Section  three  outputs  the  two  matrices  to  a  file.  A  listing  of  this  program 
can  be  found  in  Appendix  B. 

The  third  program  solves  the  eigenvalue  problem.  It  is  also  divided  into  the 
inputi  process*  output  format.  The  input  section  reads  in  the  stiffness  and  mass 
matrices  generated  from  Program  Two.  The  process  section  solves  the  eigenvalue 
problem  by  calling  the  BIGZS  routine  from  the  IHSL  library  C73.  The  output  section 
writes  the  eigenvalues  and  the  eigenvectors  to  a  file.  A  listing  of  this  program 
may  be  found  in  Appendix  C. 

The  fourth  program  generates  the  data  base  which  is  used  to  produce  a 
contour  plot  of  the  mode  shape.  It  can  do  this  as  a  result  of  solving  the 
eigenvalue  problem.  Our  admissible  functions  for  each  boundary  condition  were 
made  up  of  an  infinite  series  of  sine  functions  (each  term  the  product  of  a  sine 
function  of  x  and  a  sine  function  of  y)  multiplied  by  an  infinite  number  of 
undetermined  coefficients  (the  C  's).  When  we  solved  the  eigenvalue  problem*  the 


bottom  third  of  our  eigenvector  was  the  solution  for  the  Cmn's.  Therefore  to 
generate  the  displacement  w,  as  a  function  of  x  and  y  we  only  need  this  portion  of 
the  eigenvector  multiplied  by  its  appropriate  sine  functions  (see  Eqs.  (67),  (91), 
and  (105)).  That's  exactly  what  the  last  program  does.  It  reads  in  the  eigenvector 
for  the  desired  mode,  strips  off  the  last  third  and  then  normalizes  that  portion  of 
it.  The  program  then  generates  a  value  for  w  as  a  function  of  x  and  y  over  the 
domain  of  the  plate.  This  datafile  is  then  used  by  the  SUPBRPROC  file  on  the 
CYBER  (developed  by  Captain  Hinrichsen  and  appended  to  by  Major  Hodge)  to  plot 
the  contours  of  the  mode  shape.  A  listing  of  the  last  program  may  be  found  in 
Appendix  D. 


Several  different  characteristics  of  both  the  Qalerkin  Method  end  shear 


deformation  and  rotatory  inertia  effects  were  investigated.  To  explore  the 
Galerkin  method)  convergence  characteristics  and  comparison  to  closed-form 
solutions  were  researched.  To  study  the  effects  of  shear  deformation  and  rotatory 
inertia)  several  cases  varying  the  length  to  thickness  ratio  were  investigated.  We 
will  begin  our  discussion  by  describing  the  type  of  plates  used  in  the  subsequent 
analysis. 

Laminated  Platt  Properties 

Two  different  types  of  laminated  plates  were  used  in  this  thesis.  Goth 
consisted  of  a  graphite-epoxy  material  (AS/3501)  with  the  following  properties: 

B1  ■  21.08+06  psi. 

B2  -  1.40B+06  psi. 

v12«0.3  (116) 

Gj2  ■  0.60B+06  psi. 

P  •  0.055  lb./in.3 

where  P  is  the  mass  density  and  the  other  material  properties  have  been 
previously  defined.  One  plate  had  a  ply  layup  of  [0/903^  and  the  second  had  a 
layup  of  C±4532#.  The  thickness  of  the  plates  was  held  constant  at  one  inch  but 
the  length  and  width  dimensions  varied  (between  5  and  200  inches)  depending  on 
the  type  of  analysis  being  performed.  Tables  3.1  and  3.2  contain  the  stiffness 
elements  computed  from  Program  One  for  these  two  plate  layups. 


Graphite-epoxy  [0/90D2s 
On*  inch  thick 


enent 

Dimensional  Value 

Nondimensional  Val 

A44 

541,008 

8.385714 

A 

8 

8 

45 

**55 

548,888 

0.385714 

°11 

1,555,168 

1.11083 

D12 

35,211.3 

0.0251589 

D 

8 

0 

16 

CN 

a” 

322,770 

0.23055 

8 

0 

26 

D  . 

50.880 

0.0357143 
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Units  tor  th*  dimensional  terms  are  lb./in. 

Units  tor  th*  dimensional  term*  are  in.  -  lbs. 


Table  3.1  Stittness  Elements  tor  Plate  i 


Graphit*-*poxy  E±4532# 
On*  inch  thick 


lament 

Dimensional  Value 

Nondimensi onal 

A44 

540,808 

0.385714 

A.. 

0 

0 

45 

*55 

540,000 

8.385714 

D11 

537,089 

0.383635 

D12 

437,089 

0.312207 

°16 

308,099 

0.22007 

°22 

537,089 

0.383635 

°2  6 

308,099 

8.22007 

D66 

451,878 

8.32277 

dimensional  A^  t*rm*  *rt  lb,/in. 

Unit*  tor  th*  dimensional  terms  are  in.*  lb*. 


Table  3.2  Stittnes*  Bl»m*nt*  tor  Plat*  2 


In  using  the  Galerkin  Method  tor  solving  our  set  of  three  coupled  partial 
differential  equations,  we  have  to  ask  ourselves  some  questions.  First,  does  the 
technique  converge  to  a  solution  and  second,  how  does  the  solution  compare  with 
accepted  theory?  We  will  now  look  at  these  two  questions. 

A  discussion  of  the  proof  of  convergence  for  the  Galerkin  Method  may  be 

found  in  Reference  C9L  A  criteria  based  on  the  method  for  convergence  is  whether 

the  assumed  functions  form  a  complete  set  of  functions.  We  will  not  attempt  to 

prove  convergence  for  the  work  completed  herein.  We  will,  however,  show  the 

necessary  (but  not  sufficient)  condition  that  the  frequency  drops  by  smaller  and 

smaller  amounts  as  the  values  of  M  and  N  are  increased.  Tables  3.3  and  3.4  show 

2  3 

the  values  of  normalized  frequency  Ewa  (p/B2h  >3  for  the  three  boundary 
conditions  considered  for  the  two  plates  with  increasing  values  of  M  and  N.  The 
tables  are  based  on  a  length  to  thickness  ratio  (s)  of  20.  Boundary  condition  #i  is 
the  simply-supported  case.  Boundary  condition  02  is  the  clamped  case  and  boundary 
condition  03  is  the  clamped  (two  opposite  sides)  simply-supported  (two  opposite 
sides)  case.  The  author  could  not  compute  frequencies  greater  than  M  and  N  *  8 
due  to  computer  memory  limitations. 

As  can  be  seen  in  Table  3.3,  the  simply-supported  case  for  the  C0/9032# 
layup  converges  from  the  start.  This  is  to  be  expected  as  the  boundaries  are 
completely  satisfied.  All  of  the  other  cases  have  not  converged  but  display 
characteristics  that  makes  one  believe  they  will  converge  as  M  and  N  get  larger. 
That  is,  for  every  increase  in  M  and  N,  the  normalized  frequency  drops  by  a  smaller 
and  smaller  amount. 

A  plot  of  the  normalized  frequency  vs  M  and  N  for  BC  02  for  the  l±4532b 
plate  (third  mode)  is  presented  in  Figure  3.1.  This  condition  tends  to  have  one  of 


Non-dimensional  Natural  Frequencies  Cfc/VGDgg 
graphite-epoxy,  a/b»i,  s*20 


N 

1st  Mode 

3rd  Mode 

5th  Mode 

Boundary  Condition  Ml 

2 

11.758 

36.866 

4 

11.758 

36.866 

42.573 

6 

11.758 

36.866 

42.573 

8 

11.758 

36.866 

42.573 

2 

31.751 

Boundary  Condition  #2 
69.038 

4 

23.734 

49.403 

67.560 

6 

22.992 

48.569 

57.062 

8 

22.776 

48.328 

55.665 

2 

24.870 

Boundary  Condition  M3 
65.306 

4 

21 .235 

45.463 

52.438 

6 

20.814 

45.256 

51.832 

8 

20.697 

45.225 

51 .664 

Tabic  3.3  Normalized  Frequencies  -for  Plate  1 


Non-dimensional  Natural  Frequencies  C±453„ 
graphite-epoxy,  a/b"i,  s»20 


N 

1st  Mode 

3rd  Mode 

5th  Mode 

Boundary  Condition  #1 

2 

14.699 

36.164 

4 

14.418 

35.444 

57.883 

6 

14.283 

34.734 

55.882 

8 

14.205 

34.613 

54.856 

Boundary  Condition  62 

2 

31.091 

69.925 

4 

22.363 

44.553 

72.156 

6 

21.412 

42.964 

64.264 

8 

21.116 

42.449 

63.095 

2 

24.300 

Boundary  Condition  63 
66.437 

4 

18.956 

41 .549 

59.675 

6 

18.337 

40.438 

58.116 

8 

18.126 

40.861 

57.673 

Table  3.4  Normalized  Frequencies  for  Plate  2 


•the  "worst*  convergence  tendencies  of  the  group.  However,  even  it  does  show  very 
good  converging  behavior. 

These  tablet  are  by  no  meant  offered  at  proof  of  convergence,  but  they  do 
display  good  convergence  characteristics.  They  also  point  out  a  drawback  with  this 
method.  At  M  and  N  increase,  the  required  computer  memory  space  to  generate  the 
mass  and  stiffness  matrices  and  to  solve  the  eigenvalue  problem  becomes  quite 
large.  Thus,  if  one  wants  a  very  accurate  answer,  particularly  for  a  higher  mode, 
they  will  need  the  appropriate  computer  resources. 


The  next  question  that  mutt  be  addressed  it  the  "How  well  do  the  solutions 
generated  here  compare  to  commonly  accepted  theory?"  We  will  address  that 
question  by  investigating  a  very  long,  narrow  plate.  There  are  several  reasons  for 
this.  First,  because  of  the  nature  of  the  system  of  equations  for  our  plate,  we  can 
zero  out  the  rotatory  inertia  but  we  cannot  zero  out  the  shear  deformation  effect. 
When  shear  deformation  is  zeroed,  the  equations  become  very  ill-conditioned  and 
the  problem  cannot  be  solved.  Therefore,  we  must  look  to  commonly  accepted  theory 
that  includes  either  shear  deformation  or  shear  deformation  and  rotatory  inertia 
effects.  We  would  also  like  to  find  a  simply-supported  case  for  an  orthotropic 
plate  because  our  Galerkin  Method  converges  quickly  there.  Reference  1191  did 
contain  a  closed-form  solution  for  an  infinitely  long,  simply-supported, 
orthotropic  plate  including  shear  deformation  effects,  and  this  is  what  the  author 
used  to  validate  his  program. 


From  Reference  C 19 3  we  have 


0)  a  w'  Ci-<D,,m2t2)/(D..m2lC2  +KA,ca2)3  <117) 

mm  li  ll  55 


where 


W'  ■  (D, .  m^lt4/pa4)^ 

m  ll 


(tie) 
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H«re  w'  is  the  natural  frequency  calculated  from  tha  classical  theory  based  on  the 
re 

Kirchhoff  hypothesis.  If  we  use  the  properties  for  Plate  1  (Table  3.1)  and  compute 
the  fundamental  frequency  for  an  infinitely  long,  1ft"  wide  plate  we  have 


(i)'  >  1641.61  Hz 
re 

CO  «  1417.56  Hz 
m 


This  author  cannot  place  an  "infinitely”  long  plate  into  his  Galerkin  algorithm. 
However,  he  can  use  a  very  long  plate  such  as  one  10  inches  wide  and  206  inches 
long.  From  Reference  Cl  3,  we  can  compute  the  natural  frequency  without  shear 
deformation  and  rotatory  inertia  effects  from 

(d2p  «  Du(«14/*4)  +  2(D12-*-2D^)/(o:2/a2b2)  +  D22«a34/b4)  (120) 


where 


otj  ■  nr* 


a2  ■  niC 


2  2  A 
ag  *  m  n  1C 


for  all  m  and  n  (m  and  n  determine  the  mode). 


Again  using  the  properties  from  Table  3.1  for  our  10"  by  260"  plate  we  have 


(0  «  1,641.26  Hz 


This  compares  to  1,641.61  Hz  for  Whitney's  infinitely  long  plate  so  we  can  conclude 
that  our  260"  plate  reasonably  approximates  an  infinitely  long  one  and  we  can 
therefore  use  Eq.  (117)  to  compare  with  our  Galerkin  program  output.  Table  3.5 


shows  the  comparison  between  the  closed  formed  solution  and  the  output  from  the 


Galerkin  Method.  At  can  be  seen»  excellent  agreement  was  obtained.  Table  3.3 
also  shows  the  effect  of  rotatory  inertia  for  this  condition. 

The  author  could  not  find  information  to  compare  with  the  other  two 
boundary  conditions. 


Mode 

C 1 ose  d-Form  Solution 
with  SD 

Galerkin  Method 
with  SD  no  R1 
M  = 

Galerkin  Method 
with  SD  and  RI 

N  =  5 

1 

1417.56 

1417.56 

1414.32 

2 

1416.5 

1418.38 

1415.89 

3 

1420 . 1 

1419.97 

1416.60 

4 

1422.41 

1422.42 

1418.96 

5 

1425.51 

1425.60 

1422.08 

Note:  All  frequencies  are  in  Hz. 


Table  3.5  Galerkin  Method  and  Closed-Form  Solution 
Comparison  Plate  #1 


Shear  Deformation  and  Rotatory  Inertia  Effects 


The  second  area  of  investigation  was  to  determine  the  effect  of  shear 
deformation  and  rotatory  inertia  for  the  boundary  conditions  considered.  This  was 
accomplished  by  varying  the  length  to  thickness  ratio  and  comparing  those  results 
to  classical  laminated  plate  theory.  Only  square  plates  were  considered. 

For  the  [4/9432s  plate,  simply-supported,  the  length  to  thickness  ratio  was 
varied  from  5  to  56.  The  results  are  plotted  in  Figures  3.2  and  3.3  as  normalized 
frequency  vs  thickness  ratio  for  the  first  and  fifth  bending  modes.  It  can  be  seen 
that  shear  deformation  effects  can  be  significant  (up  to  33%  lower  for  the  first 
mode  and  52%  lower  for  the  fifth  mode  for  a  length  to  thickness  ratio  of  5)  while 
rotatory  inertia  effects  account  for  only  a  2%  lower  frequency  at  the  worst  point. 
Notice  that  as  the  length  to  thickness  ratio  approaches  56.  we  asymptotically 
approach  the  classical  laminated  solution. 

These  same  trends  are  seen  for  the  [±45^  simply-supported  plate  also. 
Figures  3.4  and  3.5  plot  the  normali'ed  frequency  vs  length  to  thickness  ratio  for 
the  first  and  fifth  modes  for  this  plate.  Although  the  frequencies  are  slightly 
higher,  the  same  behavior  occurs  as  before.  The  shear  deformation  effect  becomes 
significant  for  length  to  thickness  ratios  of  less  than  35  and  rotatory  inertia 
effects  are  very  small.  Table  3.6  presents  this  data.  There  were  no  classical 
solutions  avaliable  io  compare  these  results  with.  but.  based  on  the  behavior  of 
the  first  plate,  we  could  estimate  a  first  mode  normalized  frequency  of  about  15 
and  a  fifth  mode  normalized  frequency  of  about  63. 

Tables  3.7  and  3.8  show  the  comparisons  for  the  second  and  third  boundary 
conditions  for  the  CG/9632g  pl*te.  They  too  show  the  same  type  of  trends  as  the 
previous  two  cases.  Shear  deformation  is  a  significant  effect  below  length  to 


CLPT  -  CLASSICAL  LAMINATED  PLATE  THEORY 

SD  -  SHEAR  DEFORMATION  EFFECTS 

SDRI  -  SHEAR  DEFORMATION  AND  ROTATORY  INERTIA  EFFECTS 


i - r 

10  20 


~i - r 

30  40 


Figure  3.2  Plot  of  Normalized  Frequency  v»  Thickneee  Ratio 


Plate  #1  Simply-Supported  Boundary  Fir*t  Mode 


SHEAR  DEF0RMA1 
SHEAR  DEFORMA1 


I 


Nondimensi onal  Frequency  (i)a  (p/E^h  ) 


s 

1st 

Mode 

5th 

Mode 

SD  no  RI 

SD  and  Rl 

SD  no  RI 

SD  and  RI 

5 

9.62 

9.57 

25.44 

25.51 

10 

12.83 

12.77 

42.04 

41 .49 

15 

13.88 

13.84 

50.79 

50.39 

20 

14.31 

14.28 

55.42 

55.08 

25 

14.53 

14.51 

58.05 

57.77 

30 

14.45 

14.43 

59.45 

59.43 

35 

14.72 

14.71 

40.47 

40.49 

40 

14.74 

14.75 

41.37 

41 .23 

50 

14.88 

14.87 

42.23 

62.  i 3 

Classical  Laminated  Plate  Frequency:  None  available 

Simply-Supported  Boundary  Condition  C±45l2S  M=N=4 

Table  3.4  Shear  Deformation  and  Rotatory  Effects  for  Plate  2 


SD  no  RI 

SD  and  RI 

SD  no  RI 

SD  and  RI 

5 

10.8? 

10.85 

25.10 

24.82 

10 

17.27 

17.21 

41 .44 

41 .01 

15 

20.8? 

20.84 

51.30 

50.86 

20 

23.03 

23.00 

57.45 

57.06 

25 

24.3? 

24.36 

61 .53 

61  .21 

30 

25.32 

25.30 

65.23 

65.00 

35 

26.01 

25.?? 

68.8? 

67.  ?3 

Classical  Laminated  Plate  Frequency  1st  Mode:  26.47  5th  Mode:  74.25 
Clamped  Boundary  Condition  18/981 M =44*6 


Table  3.7  Comparison  o-f  Shear  De-formation  and  Rotatory 
Inertia  E-ffects  to  Classical  Laminated  Plate  Theory. 


Nondimensi onal  Frequency 

Ua  <p/E 

2h  > 

s 

1st 

Mode 

5th 

Mode 

SD  no  RI 

SD  and  RI 

SD  no  RI 

SD  and  RI 

5 

9.50 

9.41 

24.26 

23.96 

10 

15.25 

15.20 

36.80 

36.28 

15 

18.76 

18.72 

45.12 

44.86 

20 

20.85 

20.81 

52.09 

51  .83 

25 

22.16 

22.13 

56.97 

56.77 

38 

23.05 

23.03 

60.48 

60.31 

35 

23.78 

23.68 

63.08 

62.94 

Classical  Laminated  Plate  Frequency  1st  Mode:  24.53  5th  Mode:  70.26 


Clamped  Simply-Supported  Boundary  Condition  [0/90]^  M=N=6 

Table  3.8  Comparison  o-f  Shear  De-formation  and  Rotatory 
Inertia  E-f-fects  to  Classical  Laminated  Plate  Theory. 


thickness  ratio*  of  3®  and  rotatory  inertia  has  little  effect  on  these  bending 
modes.  These  were  not  plotted  because  a  value  of  MaN=6  was  used  to  generate  the 
data.  This  was  done  to  conserve  time  and  computer  resources.  However,  the 
frequencies  did  not  have  enough  terms  to  converge  and  at  sa50  they  overshot  the 
classical  solution  by  about  3%.  For  the  purpose  of  showing  trends  between 
rotatory  inertia  and  shear  deformation  the  author  feels  this  approach  was 
Justified. 


At  stated  in  the  ditcuttion  on  the  last  computer  programi  when  wt  tolvt 
the  eigenvalue  problem  we  alto  get  back  the  valuet  of  the  undetermined 
coefficient!  to  determine  the  mode  thape  for  a  particular  frequency  (thit  include! 
the  tortional  modet  alto).  Figures  3.6  thru  3.11  are  contour  plots  of  the  first 
modet  for  both  platet  for  ell  three  boundary  condition!.  A  length  to  thickness 
ratio  of  10  was  used  when  generating  these  plots. 


IV.  Conclusions 


Based  on  the  analysis  presented  in  this  thesis,  the  following  conclusions 

arc  presented.  They  are  organized  in  terms  of  some  general  comments  on  the 

Galerkin  Methodi  some  specific  comments  on  the  three  boundary  conditions)  and 

then  some  general  comments  comparing  the  boundaries. 

Galerkin  Method  Comments 

1.  The  Galerkin  Method  is  a  valid  approach  to  solving  the  plate  equations  of  motion 
and  yields  excellent  results  for  the  simply-supported  boundary  condition. 

2.  The  convergence  tendencies  of  the  solution  are  easily  checked  with  this  method 
by  generating  more  terms,  and  therefore  more  equations  (provided  the  assumed 
functions  are  a  complete  set  of  functions). 

3.  More  terms  were  required  to  achieve  a  converged  solution  for  problems 
considered  in  this  study  whose  assumed  functions  did  not  identically  satisfy  the 
natural  boundary  conditions  than  those  whose  assumed  functions  did  identically 
satisfy  them. 

4.  For  a  problem  that  requires  a  large  number  of  terms  to  reach  a  converged 
solution,  the  eigenvalue  problem  becomes  quite  large.  Because  of  this,  the 
author  could  not  use  more  than  8  terms  in  the  Galerkin  equations.  This  effected 
the  quality  of  results  for  the  clamped  and  clamped  simply-supported  boundary 
conditions. 

Simply-Supported  Boundary  Comments 

i.  The  natural  frequency  for  a  square  orthotropic  plate,  simply-supported,  for  the 
first  mode  is  !•%  lower  than  the  Classical  Laminated  Plate  Theory  (CLPT) 


•frequency  at  a  length  to  thickness  ratio  (represented  by  s)  of  11  And  drops  to 
33%  lower  for  An  s  of  5.  For  the  fifth  mode,  the  frequency  is  10%  lower  At  An  s 
of  23  And  drops  to  52%  lower  for  An  s  of  5. 

2.  For  a  cross-ply  laminate  with  the  seme  conditions  as  Above,  the  seme  type  of 
sheer  deformation  effects  were  seen.  No  CLPT  solution  was  found  in  the  current 
literature,  therefore  no  comparison  could  be  made. 

3.  For  both  types  of  plates  above,  the  analysis  of  rotatory  inertia  effects  showed 
a  lowering  of  the  natural  frequency  from  the  shear  deformation  frequency  by  one 
to  two  percent  for  s's  below  25  for  the  first  five  modes. 

Clamped  Boundary  Comments 

1.  The  natural  frequency  for  a  square  orthotropic  plate,  clamped  on  all  sides,  for 
the  first  mode,  is  10%  lower  than  the  CLPT  frequency  for  an  s  of  21  and  drops  to 
59%  lower  for  an  s  of  5.  For  the  fifth  mode,  the  frequency  is  10%  lower  at  an  s 
of  33  and  drops  to  66%  lower  for  an  s  of  5. 

2.  For  the  same  case  above,  the  value  of  the  asymptotic  limit  for  the  natural 
frequency  was  3%  higher  than  the  CLPT  solution.  This  was  due  to  the  fact  that 
an  M  and  N  of  six  was  used  to  generate  the  data.  To  achieve  closer  results,  a 
larger  value  of  H  and  N  would  have  to  be  used. 


3.  The  analysis  of  rotatory  inertia  effects  showed  a  lowering  of  the  natural 
frequency  from  the  shear  deformation  frequency  to  be  at  mast  1%,  occuring  at  an 


1.  The  natural  -frequency  for  a  square  orthotropic  plate,  clamped  on  two  opposite 
sides,  simply-supported  on  two  opposite  sides,  for  the  first  mode,  is  10%  lower 
than  the  CLPT  frequency  for  an  s  of  25  and  drops  to  61%  lower  for  an  s  of  5. 
For  the  fifth  mode,  the  frequency  is  10%  lower  at  an  s  of  35  and  drops  to  65% 
lower  for  an  s  of  5. 

2.  For  the  same  case  above,  the  value  of  the  asymptotic  limit  for  the  natural 
frequency  was  1.5%  higher  than  the  CLPT  solution.  This  is  again  due  to  the  fact 
that  a  value  of  M  and  N  equal  to  six  was  used  to  generate  the  data. 

3.  The  analysis  of  rotatory  inertia  effects  showed  a  lowering  of  the  natural 
frequency  from  the  shear  deformation  frequency  to  be  at  most  1%,  occuring  at  an 
s  of  5. 

General  Comments 

1.  The  effects  of  shear  deformation  are  more  significant  for  the  two  clamped 
boundary  conditions  than  for  the  simply-supported  boundary. 

2.  The  effects  of  shear  deformation  increase  with  increasing  mode  for  all  three 
boundary  conditions. 


3.  Analysis  shows  that  rotatory  inertia  has  very  little  effect  for  all  three  modes. 
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•  9  REM  THIS  PROGRAM  COMPUTES  THE  EXT EN SIONAL  AND  BENDING 
,  REM  STIFFNESS  ELEMENTS  FOR  A  SYMMETRIC  LAMINATE  BUILD 
30  REM  UP,  GIVEN  LAMINA  PROPERTIES. 

40  REM 
50  A 1 =0 
60  A 2=0 
70  A 3=0 
80  D1=0 
90  D2=0 
100  D3=0 
110  D4=0 
120  D5=0 
1 30  D6=0 
140  PS=0 

150  INPUT* DO  YOU  WISH  TO  ADD  ANOTHER  LAYER* ;N$ 

160  IF  N*=*N*  THEN  GOTO  860 
170  REM 

180  REM  THIS  SECTION  GETS  THE  LAMINA  DATA. 

190  REM 

200  INPUT "ORIENTATION  ANGLE  " ;TD 
210  TH=TD*3. 1415927/180 
220  INPUT* 2K  DIMENSION  *;ZK 
230  INPUT* ZK-1  DIMENSION  ";Z1 

240  INPUT"ARE  THE  REST  OF  THE  LAMINA  DATA  THE  SAME  AS  THE  LAST  TIME  * ;M* 

250  IF  M*=*Y“  THEN  GOTO  440 
260  INPUT* El  * ;E1 
(1  0  INPUT* E2  *  ;E2 
280  INPUT* VI 2  * ;V1 
290  V2=V1 *<  E2/E1 ) 

300  INPUT* G1 2  * ;G1 
310  G3=G1 
320  G2=0.8*G1 

330  INPUT’MASS  DENSITY  * ;RHO 
340  REM 

350  REM  COMPUTE  THE  A  AND  D  ELEMENTS  FOR  THIS  LAYER. 

360  REM 

370  Ql=El/< 1 -< VI *V2) ) 

380  Q2=(Vl*E2>/< 1-<V1*V2)> 

390  Q3=E2/< 1 -< VI *V2) ) 

400  Q4=G2 
410  Q5=G3 
420  Q6=G1 
435  REM 

440  REM  COMPUTE  THE  QBARS. 

445  REM 

450  B1=Q1*(C0S<TH))[4  +  2*< Q2+2*Q6) *< SIN<TH> > [ 2  *(C0S<TH)>[2  +  Q3*( SIN<TH) > C 4 
168  B2=<  Q1 +Q3-4*Q6) *<  SIN(TH) ) [ 2  «  <C0S<TH>>[2  +  Q2*( <  SIN<TH) >14  *  < COS(TH) ) C 4) 
470  B3=Q1*<SIN<TH) ) [4  ♦  2*<  Q2+2*Q6) *<  SIN<TH> ) C 2  *  <C0S<TH)H2  +  G3*< COS<TH) ) [ 4 
480  B4=<Q1-Q2-2*Q6>#SIN<TH)*<C0S<TH> ) C3  +  <  Q2-Q3+2*Q6>  *<  SIN(TH) > 1 3*C0S<TH> 

490  B5=<Q1-Q2-2«Q6>*<SIN<TH> ) C3*C0S<TH)  ♦  < Q2-Q3+2*Q6> *SIN<TH> *< C0S<T:O  > [ 3 
508  B6=<  Q1 ♦Q3-2*Q2-2*Q6) *<  SIN<TH>  >  t  2*<  COS(TH) >12  +  Q6*< ( SIN<TH> > 1 4  ♦  < COS<TH> > C 4 
). 

J  B7=Q4»<  COS<TH)  )  t  2  ♦  Q5*<SIN<TH>>[2 
520  B8=<Q4- Q5)*C0S<TH)*SIN<TH> 


530  B9=Q5*(C0S<TH) > [2  +  Q4«< SINCTH) ) [ 2 
550  REM 

560  REM  COMPUTE  THE  A  ELEMENTS. 

570  REM 

580  A4=B7*<  ZK-Z 1 ) 

590  A5=B8*<  ZK-Z  1  ) 

600  A6=B9*<  ZK-Z  1 ) 

610  P=*RHO*<  ZK-Z  1  ) 

620  REM 

630  REM  COMPUTE  THE  D  ELEMENTS. 

640  REM 

659  DZ=<  ZKC  3-Z 1 [ 3) 

660  Fl=< 1/3)*B1*DZ 
670  F2=< 1/3) *B2*DZ 
680  F3=< 1/3)*B4*DZ 
690  F4=< 1/3>*B3*DZ 
700  F5=< 1/3>*B5*DZ 
710  F6=< 1/3) *B6*DZ 
720  REM 

730  REM  SUM  THE  A/ S  AND  D'S  WITH  THE  PREVIOUS  LAYERS 

740  REM 

750  A1=A1+A4 

760  A2=A2+A5 

770  A3=A3+A 6 

780  D1=D1+F1 

790  D2=D2+F2 

800  D3=D3+F3 

810  04— D4+F4 

820  D5=D5+F5 

830  D63CD6+F6 

840  PS=PS+P 

850  GOTO  150 

854  REM 

355  REM  PRINT  OUT  THE  A  AND  D  ELEMENTS. 

856  REM 

860  LPRINT" A44  -  * ;A1 
870  LPRINT* A45  =  “ ;A2 
880  LPRINT-A55  =  * ;A3 
890  LPRINT* D1 1  =  * ;D1 
900  LPRINT“D12  ■  ’ jD2 
910  LPRINT* D1 6  *  * ;D3 
920  LPRINT* D22  =  * |D4 
930  LPRINT* D26  =  *;D5 
940  LPRINT* D66  =  *;D6 
950  LPRINT* P  =  * ;PS 

960  REM  NOW  COMPUTE  THE  NORMALIZED  STIFFNESSES. 

970  INPUT* INPUT  THE  PLATE  THICKNESS* }H 
980  A1»*A1/<E2*H> 

990  A 2= A 2/ < E2*H) 

1000  A3=A3/<E2*H> 

1010  Dl=Dl/<  E2*H£  3) 

1020  D2=D2/ ( E2»H£  3) 

1030  D3=D3/<  E2*Ht  3) 

1040  D4=D4/<  E2*H£  3) 

101 


1058 

I860 

1070 

1080 

1090 

1100 

1110 

1128 

1130 

1140 

1150 

1160 

1170 

1180 


D5=D5/< 

D6=D6/< 

LPRINT 

LPRINT* 

LPRINT* 

LPRINT* 

LPRINT* 

LPRINT* 

LPRINT* 

LPRINT* 

LPRINT* 

LPRINT* 

LPRINT 

END 


E2*Hl 3) 
E2*H[ 3> 


*  ;A1 

*  ;A2 
"  ;A3 

*  ;  D 1 

*  ;D2 

*  ;D3 

*  ;D4 

*  ;D5 
“  ;D6 
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r>  n  n  o  o  r> 


PROGR AN  8CTW0 <  I MPUT »OUT PU T »T APE 5  = INPUT « TAP 1 6  =  OUTPUT ) 
REAL  Fl»F2fF3*F4,F5,F6»F7»F8,F9 ,F10 »K tK4,K5«Pl 
INTEGER  MfNtP*Q«ITESTl.ITEST2»NNAX*NHAX 

REAL  AlC10*10>»A2C10tlO)tA3<10tlO)#BlC10*10)»82<10*10> 
REAL  BJ(10»lG)*Cl(10»10)»C2(10*10)»C3(10flO> 

REAL  IN(108»10e)»FCC108>*STC108  »108 )*PAC108»108) 

INPUT  DATA  FOR  THE  CONPCSITE  PLATE. 

MHAX=b 
NMAX=HNAX 
PI=3. 1415927 
R  =  1.0 
A=10.0 
B  =  A 
H  =  1  .0 
S  =  A  7  H 

Pl=1.4245£-C4 

K=5./6. 

K4=1.0 
K  5= 1 . 0 

A44=3.«571F-01 

A55=3.8571E-01 

Dll=l  .11083 

D12-0 .9251509 

016  =  0  .  0 

022=0.23055 

026=0.0 

OO,  6  =  0.  0357143 

ET=1. 45*06 

DETERNIM  THE  VALUE  OF  THE  INTEGRATED  TERNS. 


1  =  1 
J  =  1 

00  190  P=1*NHAX 
DO  109  G= 1 »  NHAX 
DO  50  N= 1 » HHAX 
DO  50  N=1,KNAX 
ITE$Tl=  HOD  <N*P  »2  ) 

I  TE  ST2=  HC0(N-*Q.2> 

IF  CN.EQ.P)  THEN 
F 1  =  C .5 
F  3  =  0 • 0 
F  4  =  C  •  0 
ELSE 

F 1  =  0  •  0 

IF  < ITEST1.EQ.0)  THEN 
F  2  =  0 • 0 
F  4  =  0 • 0 
ELSE 

F3=2*N 

F4=2*P 

ENOIF 

ENCIF 

IF  (N.EQ.Q)  THEN 


it 


r2  =  3.0 

IF  (ITEST2.EC.O 

F  =  =  3  •  0 
F  f .  =  0  .C 
ELSE 

F'5  =  2«N 
Fb=2*G 
ENOIF 
END  IF 

F  7= (H*F  )- <P«PJ 

F8=  -F7 

IF  (F7.FQ.0.)  THEN 
F7  =  l. 

F  9=  1  • 

ENDIF 

F9  = (N*N)- <Q*Q) 
F10=-FS 

IF  (F9.EQ.0.)  THEN 

F9=l  . 

F 10  =  1  • 

ENDIF 


THEN 


COMPUTE 


STIFFNESS  MATRIX. 


A1  (H#N>=PI»PI*F1*F2M-M*H*D11-MMR*R*''6E>-K*S*S*A?5/(PI  *PI>  ) 
1*FA*F6* (2«N«N*R*D16 )/(Ffi*F10) 

B1  <H,N)-PI*PI«F1*F2»  <  -M*H«D16-h*MR*R*C26> 
1*FA*F6.<H*N*R*<D12-»C66  )  )/«F8*FlC) 

Cl CM,N)=-FA*F2* {M*K«S*A55>/F8 
A2(f',N)=PI*PI*Fl*F2*<-M*H*D16-MMR*R«D26> 

1*F4*F6*  (M-N*R*<012-»C66J)/CF8*F1C> 
B2<M,N>=PI«PI*F1*F2M-M*H*066~MN<*R*R*D22-K«S*S*A4  4/(PI*PI> ) 
l>F4*Ff •C2*M*N*R*D26/<F8*F10) > 
C?(MtN)-”Fl«F6*<N*R*K*S*A4A/F10) 

A3(MtN)=F4*F2*CM*K*S*A55)/F8 

R3tM,N>  =  Fl«F6*<N*R*l'*S«A4A)/F10 

C3 <M  .Nl -P I*PI*F1*F2 • (-M*M«K*A55-N*N*R*R*K* A 44) 

ST<  I.d)=-A1(M,N) 

ST  <  I  »J*MMAX*NMAX)=-E1(M»NI 
ST( I »d«2*MMAX*NMAX»  =-Cl <HtN> 

ST<  I*MMAX*NMAXt J>=-A2<M*N) 

ST< I*MMAX*NMAX*J»MHAX*NMAX>  =  -B2  CM|N> 

ST<  I-*MHAX*NMAX»d*2*MHAX*NHAX>  =  -C2 (H»N) 
ST<I^2»MMAX*NMAX»J>=-A3TMtN> 

3T(  I*?«MHAX*NMAX»  J*MMAX*NMAX  >  =  -E3<MtN> 

STC  I*2«MNAX*NMAXf  J»2*MMAX*NMAX>=-C3<N*M 

J  =  J*l 

CONTINUE 

I=I«1 

J  =  1 

CONTINUE 


COMPUT 


HASS  MATRIX 


DO  20  1  =  1  «  <3*MHAX*NNAX  ) 
DO  20  «J  =  1  *  (3*HHAX*NMAX  I 
IF  (I.EQ.J)  THEN 
IM( I tJ  1  =  1.0 

ELSE 

IH(I »  J )  =  0 • 0 


v\ 


•«  L  -  A 

V‘-V-' 


*  k  *  »  *  h  '  m  •  t.  “  >  ‘  .  ‘  "  W  *  .  -  .  *  •  *  .  •  «  *  J.  \ »  *  W  '  >  *  ■ 


■  v.VV 


END  IF 

20  CONTINUE 

DO  30  1  =  1  t  (2*HHAX*NMAX  )  i 

MC<I)=K‘;*<P1*A*A)/(A8.C*H*ET) 

3  C  CONTINUE 

00  A  3  I  =  (  2*NMAX*NHA  X  *1  )  <  (3*HMA  X  »h  HA  X  ) 
MC(I)=M*(<R*,*A)*<R**A)*P1)/(A»C*S*S*ET*(H**3)) 

AC  CONTINUE 

DO  £0  1  =  1  ,  (3*MMAX*NMAX>  J 

00  60  J-l » <3*MMAX*NMAX > 

MA(I,J)  =  IMIt.J>*MC(0> 

6  C  CONTINUE 

WRITE  THE  STIFFNESS  AND  HASS  MATRICES  TO  A  DATA  FILE. 

w* 

DO  151  I=1,3*MMAX*NMAX  *. 

00  1  a  1  J=1 ,3*MMAX*NFAX 
WR  ITE(fc«150l  STCI »d  >»MA<ItJ) 

150  F0RMAT(2E1 3.5) 

151  CONTINUE 
STOP 

ENG  c 


\  W 

i 


Equations  for  Boundary  Condition  #1 


COMPUTE  THE  STIFFNESS  MATRIX. 


A1  (li»N)-P  I*FI*F1*F2*<-M*M*011-N*N*R*R*066-K*S  *S  *A  55/C  PI  *P  I  >  > 
1-F3«F6*(2*R*M*N*D16  )/<F7*F10) 

2*F6* (N*R*D16*  C1-COS<M*PI>*COS(P*PI)  )  >/F10 
81 ( WfN»=PI*PI*Fl*F2*{-M*N*R*  CD15*D66) I 

l-F3*F6*<M*M*O16+N*N*R*R*D261/(FT*F10) 
2*F6*(M*D16*(1“C0S<M*PI)*C0SCP«PI) )) /F10 


C1(M«N)=-PI*F1*F2*M*K*S*A55 
A2  <H»N)=PI*PI*F1*F2*(-H*N*RMD12*D66>  ) 
1-F4*F5«(M*M*D16*N*N*R*R*026>/{F8*F9> 
2.fa*(N*R*026*(1-COS(N*PI)*COS(Q*PI) >>/F8 
B2(H»N)=PI*PI*F1*F2*C*,H*M*D66-N*N*R  *R*022-K*S  •S*A44/<PI*PI>> 


l-F4*Fb*(2*M*N*R#D26/<F8*F9) 1 
?.F4*M*026*( l-COS< N*FI)*C0S<Q*PI)>/F8 
C2(M  *N)--PI*F1*F2*< MR*K*S*A44I 
A3( H,N)=-PI*F1*F2*M*K*S*A55 
R3(M»N)=-PI*F1*F2*N*R*K*S*A44 

C3(M,N>=PI*PI»F1*F2*(-M*M*K*A55-N*N*R*R*K*A44) 


ST(I»J>=-A1<H.N) 

ST< I » J*MMAX*NHAX>=-BltN»N> 

ST  < I ,J*2*MMAX*NMAX>  =-Cl(M»N) 

ST< I^MMAX*NMAX* J>=-A2CMtN) 

ST ( I*MMAX*NMAX» J*HM AX*NMA X > =-B 2 < M »N > 
ST<I*MMAX*NMAX*J*2*MMAX*NHAX»=-C2(M|N> 
ST<  I«’2*HMAX*NMAX*  J)  =-A3(MfN) 
STCI»2»MMAX*NMAX»J*MMAX*NMAX)=“B3<MtN) 
ST(I»2*HMAX*NMAX*  J*2*HMAX*NHAX>=-C3(M,N> 
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• 


n  n  r> 


.-v  - 


► 

I 


Equations  for  Boundary  Condition  ??3 


COMPUTE  THE  STIFFMESS  MATRIX. 

A 1  O*#N)=PI*PI«Fl»F2«<-M*M011-MMR*R*Cb6-K*5*$*A55/CPI*P 
1*F4*F6*C2*R*M*N«016  >/<Ffi*F10> 

B1 CM»N)~-F1*F£*PI*<  MM*Cl€*h*N*R«R*G26  )/FI0 
l“FA*F2*PI*<M*N*R*«012-»C£6n/F8 
Cl CH,N)=-F4»F2*H*K*S*A5S/F8 
A2CM*N>=:-F1*F5*PI*<  F*M *01«*K*N* R* R* 026 > /F9 
l*F4*F2*PlMM*fv*ft*<D12*066>WF8 
2*Fl*PI*N*R*02  6*a-CCS(N«PI>*CC5(a*PI>  ) 
B?CM,N»=PI-PI«F1*F2«<-M*K*066-K*I»*R*R*D22-K*S*S*A44/<PI*P 
1-'F4*F5*(2*M*N*R*026/(F8*F9) > 

2«'F4*H*Q2‘i«  ( l-COS(N*FI  )*CQS(Q*P1»)/F8 
C?<MtN>  =  -PI*Fl*F2*<A*R*K*S*AM) 

AT  CM  tN)=FA»F2*M*K*S'*A5  5/F8 
B3(M»N>r~PT*Fl*F2*K*R*X*S*AM 

C3CF»N>=PI*PI*F1*F2«(-MM«K* A55 -K *N *R« R «K* A 44  ) 
ST(I«0>=-A1CM,N> 

ST(I»J-*MMAX*NPAX)  =  -E1<M»N) 

ST  C  I  tJ-*2*MMAX*NKAX>=-Cl  t*»IO 
STC  I*HHAX*NMAXf0>=-A2OSNI 
ST<I>HWAX*NMAX» AX* NMA X > =-B 2  <K*N> 

ST<  I*MI«AX*NMAX»*J*2*FFAX*N*AX>  =  -C2  CB»*> 

ST(I*2*MHAX*NMAXf J)=-A3CR»N> 

ST<  I*2*XMAX*NMAX»U'*MMAX*NKAX>  =  -E3(HiN> 

ST  U*2*MMAX*NMAX»  J* 2 *MK AX* W*A X > =-C3 C H » M 


Appendix  C 


PROGRAM  I  G  F I  N  (  INPUT tCUTPUT  *  T  A  F  r  S  =  I  NP  0  T  *TAFLF> -  OUTPUT) 
•internal  :ir.7s 

R 1  u  M  -t  >  ♦  »i  C  7t»  >  »  i  (  1  2  »  1  :  )  »  C  (  1  2  )  »  W  K  (  1  5  8  ) 
o  A L  ST(  »3*\irt) tVA(9e«3H) 
i n t  _ s  r •;  i  ,  j  ,r  ,l  .mm  a  >  ,n  ,  ijcb  » iz  *  ier 
C 

M  M  A  X  -  ? 

C 

C  L  ■ '  A  J  IS  T  H :  MAGS  AN  G  STIFFNESS  MATRICES. 

C 

1)0  1  "  !  -  l  »  N 
T.  1  '  U  =  1  .  N 

R  .  A  :  (  '  ,  "  )  *.T  (  I  tJ)  .  MA<  I  ,J) 

2C  K  JR  •' AT  <  .  -  1  3  .r> ) 

1C  C  \  '  T  M. 

c 

C  PUT  TH  STIFFNESS  AND  MAGS  MATRICES  IN  SYMMETRIC  FORM 

C 

T  -  i 

Dl  ‘j  '  '  i  i  N 

m  7  :  E  -  1  ,  c 

A  (  T  )  -  G  T  (  R  .  L  ) 

:  t(  I  i  C  ) 


c  :  TH-  f  IG*  NV  UUr.S  AND  EIGENVECTORS. 

C 

»  y  -  f, 

i  >i : u  " 

CALL  -  !  ./3(A,P,N,I0CBtC»Z»I7*WK.lLR) 

UR  I  T  (  ,  '*  '  )  I  c.R 

A  '  F  Jr>  M  A  T  (  I  A  ) 

■j  i  ■  ;  - 1  t  n 

UR  I  Ti-  (■»-.:  )  I «0<I > 
b  C  f  a  R  *  A  T  (  •  0  <  »  »  I  2  .  •  )  =  ' »E 12* -  ) 

5  C  C  'NT! NL. 

D  V  ■'  I  -  1  .  N 
J  3  /  n  0-1  tN 

us  i  tl  ( ,3 : )  i  iJfZd  »*j  > 

Pr  FORMAT ( *7 *«I2fiy*I2**  =  »,E12.5> 

7  r  C  1 N  T  I  M  •. 

ST  CP 
L  N  C 
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«  >_  4  5  "r  f  INPUT  »  CUTPLT  »  T  AP  ■  _  =  I NP  LT  »TAP^^"OLTPoI) 
K--  M  x  ,  1 ,  ?  ,'I  ,K  (1  )  ,C  (  i. . G )  *NC  C  i  ,  ‘  )  .LARGr. 

I  '4  T  S  I  i  J  <  L  i  C  t  N 

3  i  -  ;  7 

A  -  *  •  i 

v  =  -  . : 

T  I4  I  J  cna  C  -  T  S  TP-  C  (  X  *  N  )  CCGFFICItNTS 

O’''  1  J  L  -  1  »  1  F'. 

H  .  V" ( r  .  1  -  >  K  (  L  ) 

F  X  P  A  T  (  -  :  <  .  :  ) 

C  ;  N  r  i  N  L 
!  =  ?  ■» 

D' •  *■  - ; « 

r».  '  X  r  >  ,  . 

C  ( iJ  »  )> -x  (  I  ) 
i  =  i  * : 

C  ■ 1  f.  T  I  '•  G 

f)Vj  I  TP-  C(H,K)  CCfFFICItNTo 

L  A  •'  r  :  =  i  (C  <  1  »  :  >  > 

D''  :  s  v  = :  * ' 
i  >  ■  .\  = : » • 

I  F  (  JX  :  (  C  (  *  ,  f. )  )  .  0  T  .  L  i  *  C  >  T  P  :  N 
L  *  K  j  =  a  -  •-  ( c  (  *  .  \  )  ) 

:  •  i  f 
C  T  1\  L 
7 ' 

n  •.  :  '>  •  = :  *  r 

N  G  (  P  ,  N  >  '  G  (  P  *  f  J  >  /  t  A  «  G  r 
c \  ,■ :  i 

r h i '  s  cri  tN  c^DuT;:  tp:  platf  c ispl acfpent  « 

j  '  « '  '  = 

x  =  <  ♦ „ .  i 

D ‘  o-l  »  2  2 
v  r  /  *  i  .  * 


[i  i  . 1 

d  i  '■  n  - :  » - 

7  -■’•»(  N  C  (w»N)*SIMP*ri*X  )*SIN(N-*FI*Y  )) 

c:\tim 

vjr  n-;i  •  » r :  >  x  »y  i? 

C  ■.  N  T  I  M 
v  =  -  . : 

C  JM!  M  • 
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